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Abstract 

We compute the master integrals containing 2 and 3 massive propagators en- 
Q^' tering the planar amplitudes of the 2-loop electroweak form factor. The masses 



Oh! 



of the W, Z and Higgs bosons are assumed to be degenerate. This work is a 



D , continuation of our previous evaluation of master integrals containing at most 

' 1 massive propagator. The 1/e poles and the finite parts are computed ex- 

actly in terms of a new class of 1-dimensional harmonic polylogarithms of the 
^ . variable x, with e = 2 — D/2 and D the space-time dimension. Since thresh- 

H ' olds and pseudothresholds in s = ib4m^ do appear in addition to the old ones 

in s = 0, ibm , an extension of the basis function set involving complex con- 
stants and radicals is introduced, together with a set of recursion equations 
to reduce integrals with semi-integer powers. It is shown that the basic prop- 
erties of the harmonic polylogarithms are maintained by the generalization. 
We derive small- momentum expansions |s| <C m^ of all the 6-denominator 
amplitudes. We also present large momentum expansions |s| ^ of all 
the 6-denominator amplitudes which can be represented in terms of ordinary 
harmonic polylogarithms. Comparison with previous results in the literature 
is performed finding complete agreement. 
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1 Introduction 



We present in this paper the computation of the master integrals containing 2 and 3 
massive propagators entering the planar amplitudes of the 2-loop electroweak form 
factor. This work is the natural continuation of our previous evaluation of master 
integrals containing at most 1 massive propagator. The model process we consider 
is: 

/(Pi) + /(P2)-X(g), (1) 

where // is an on-shell massless fermion pair, Pi = P2 = ^, and X is a singlet under 
the electroweak gauge group SU(2)l x f/(l)y. 

At the 2-loop level, the annihilation in Eq. (Q) involves the emission of 2 virtual 
bosons among 7, W , Z and if's. The cases of emission of {€) 2 photons and {ii) 
1 photon and 1 massive boson, have already been treated in our previous work ^J. 
Since the above amplitudes have only thresholds in s = 0, m? and pseudothreshold 
in s = — m^, we succeeded to express them in terms of 1-dimensional harmonic 
polylogarithms j2| El with maximum weight up to 4 included. 

In this paper we consider the amplitudes with 2 or 3 bosons of mass m exchanged: 
the masses of W , Z and Higgs bosons are assumed to be degenerate, mw ^ niz ~ 
mn ~ "m..^ In general, these amplitudes have thresholds in s = 0,m^,4m^ and 
pseudothresholds in s = — m^, —Am?. As a consequence of this more complicated 
structure, the basis function set of the harmonic polylogarithms considered in |2j 
is no more sufficient and a generalization is presented. New basis functions such 
as 1/(4 ±x), l/iyx(T±~r), etc. are introduced and recursion relations to reduce 
integrals with semi-integer powers coming from the evaluation of the master inte- 
grals are derived. The basic properties of the harmonic polylogarithms, such as 
the uniqueness of representation as repeated integration, the algebra structure, the 
closure under the inverse transformation x — etc. are all maintained. 

We use dimensional regularization |^ to regulate both ultraviolet and infrared 
divergences, which then appear as poles in 1/e, with e = 2 — D /2 and D the space- 
time dimension. The ultraviolet poles are related to coupling constant renormal- 
ization and are subtracted with ordinary renormalization prescriptions, such as for 
example the MS scheme. The infrared poles are not physical, and in the physical 
cross-section, are canceled by the corresponding ones appearing in the real photon 
emission contributions or are factorized into QED structure functions. The use of 
dimensional regularization as double regulator for both infrared and ultraviolet di- 
vergences makes a priori impossible to trace the dynamical origin of the 1/e poles. 
This is however possible by using power counting estimates for the ultraviolet and 

^Small corrections of order g'^{'m\ — rriy^)'^ can be included by means of expansions of the 
denominators of the form 

The amphtudes with powers of the W denominators can be reduced to the master integrals by 
means of the integration-by-parts identities (see Section 2). 
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the infrared regions together with some physical intuition (see sec. (0}). In the ul- 
traviolet region one can usually neglect masses and external momenta while in the 
infrared region massive lines can be shrunk to a point. Furthermore leading 1/e 
poles often originate from ordered regions such as kf <^ k^, with ki and k2 the loop 
momenta. This (rather qualitative) analysis offers some checks of the results. 
The paper is organized as follows. 

In Section 121 we outline the strategy for the exact analytical evaluation of the 2- 
loop electro-weak form factor. The presentation is rather sketchy and we refer to our 
previous work for a detailed discussion of the various steps of the computation. 

In Section El we discuss the properties of the harmonic polylogarithms. An 
overview of the ordinary harmonic polylogarithms is presented, in order to intro- 
duce later the generalization necessary to represent the master integrals. As already 
anticipated, the basic tool is a set of recursion equations to reduce the integrals 
coming from the evaluation of the master integrals to a unique form, specified by 
the choice of the basis functions. 

Section jH contains the main results of our work. We present exact analytical 
results for the 21 non-trivial master integrals containing 2 and 3 massive denomi- 
nators. All the master integrals are represented in terms of generalized harmonic 
p olylogar it hms . 

In Section El we give the results for the reducible 6-denominator diagrams. The 
latter are 3 amplitudes of vertex-insertion type (see later), which are interesting by 
themselves or for reference use. 

In Section ini the small momentum expansion |s| <^ of all the 6-denominator 
amplitudes is presented. We compare our results with those present in the literature 
finding complete agreement. 

In Section [71 we compute the large momentum expansion | s | ^ of all the 
6-denominator amplitudes which can be represented in terms of ordinary harmonic 
polylogarithms. We could not give the large momentum expansion of the amplitudes 
containing generalized harmonic polylogarithms H{td; x), because the expansion of 
the latter for |x| ^ 1 is still an open problem. 

In Section IHl we draw our conclusions and we discuss the open problems in the 
computation of the master integrals. 

In order to make the paper as clear as possible, we added 4 appendices. 

In Appendix^ we report the 1-loop master integrals with 2 massive propagators 
which occur in our computation. They are written in terms of generalized harmonic 
polylogarithms, like the 2-loop amplitudes. 

In Appendix^ we present the 4 factorized 2-loop master integrals, which are the 
product of one-loop master integrals. Their derivation is rather trivial and we give 
them for completeness. 

In Appendix (Owe compute some interesting 2-loop amplitudes containing up to 
5 denominators included which can be reduced to simpler topologies by means of 
the integration-by-parts identities. 

Finally, in Appendix O we give the reducible 6-denominator scalar amplitudes 
of self-energy insertion type. 
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Figure 1: Vertex-correction diagrams with 2 and 3 massive propagators. The 
topologies related to these diagrams are real 6-denominator topologies (see text). 
The graphical conventions are the same as in our previous paper p]. 
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Figure 2: Self-energy correction diagrams with 2 and 3 massive propagators. The 
topologies related to these diagrams are 5 and 4 denominator topologies (see text). 
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2 The computation of the form factor 



The planar diagrams containing 2 and 3 massive propagators entering the 2-loop 
electroweak form factor are plotted in Figs. ^ and El We have omitted 1-particle- 
reducible contributions, related to field and mass renormalization of the external 
particles /, / and X. The amplitudes can be obtained, for example, by means of 
a Dyson-Scwhinger expansion of the ffX vertex function Let us look at the 
general structure of the diagrams: 

• 1-loop: one has a triangle diagram with a vector particle, i.e. a 7, a or a 

exchanged between the fermions; 

• 2-loops: the triangle diagram has its bare propagators and vertices replaced 
by the corresponding one- loop quantities^. 

According to the kind of correction, the 2-loop diagrams are then naturally classified 
as: 

1. Vertex-correction amplitudes (see Fig. Q). 

They all involve 6 different denominators, the maximal number of denomina- 
tors for 2-loop 3-point functions. According to the kind of vertex corrected, 
there are the following cases: 

(a) Ladder diagrams. 

The "hard" ffX vertex, where annihilation occurs, is corrected. These 
are the diagrams (a) and (b) in Fig. ^ 

• diagram (a) represents the exchange in the t channel of a pair of W/Z 
particles between the fermions and is a "non-singlet" contribution to 
the process; 

• diagram (6) represents the conversion of the initial fermion pair into 
a boson pair which convert back into a fermion pair; fermions are 
emitted in the t-channel. This diagram is a non-singlet contribution 
which actually cannot be obtained as correction of the basic one-loop 
amplitude. 

(b) Vertex-insertions. 

One of the 2 vertices involving the external fermions is corrected. There 
are 5 of such diagrams: 

• diagram (c) containing an abelian correction with a W or a. Z to the 
basic one-loop amplitude; 

• diagram (d) and (e) representing respectively the emission by an in- 
termediate boson of a photon internally or externally to the triangle; 

• diagram (/) containing the conversion of a photon into an interme- 
diate boson pair, i.e. the process 7* W*W*] 

^We neglect in this paper the crossed ladder topology, which is a non-planar one. 
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• diagram (g) containing the splitting Z* — > this diagram is 

the only one having 3 massive propagators. 

The vertex-insertion diagrams are on the same footing as the ladders but 
are conventionally named in a different way. 

2. Self-energy correction amplitudes (see Fig. |21). 

There are different cases, according to the kind of line corrected and the type 
of interaction. Self-energy-corrections may indeed occur on: 

fermion line. 

This is the case of diagram (b) in Fig.|21 with a bubble insertion. This am- 
plitude has 5 different denominators, the one representing the corrected 
line being squared; 

(b) a boson line. 

This case is more complicated than the previous one because one can 
have mixing and quartic interactions. Concerning the first point, one can 
have: 

i. a "diagonal" correction. 

By this we mean a correction to a 7, ly or Z propagator, involving 
a virtual fermion or boson pair. These are the cases (6), (c), (e), (/) 
and (h) in Fig. |21 

ii. 7 — Z° mixing. 

By this we mean a contribution to the 7 — Z° propagator. In this 
case one has 2 propagators with the same momenta but with different 
masses, which can be disentangled with partial fractioning pp. These 
are the cases (d) and (g) in Fig. |21 

Let us now consider the second point. Because of the presence of both 
cubic and quartic interactions, one can have: 

i. bubble insertion. 

These are the cases (a) — (/). There are 5 different denominators. In 
the "diagonal" case, the denominator representing the corrected line 
is squared. In the 7 — Z° mixing case, the amplitude is a superpo- 
sition of two 5-denominator amplitudes, one containing the photon 
propagator and the other the Z'^ propagator; 

ii. tadpole insertion. 

These are the cases (g) and {h). There are 4 different denominators. 
In the "diagonal" case, the denominator representing the corrected 
line is squared. In the 7 — Z° mixing case, the amplitude is a super- 
position of two 4-denominator amplitudes, one containing the photon 
propagator and the other the Z° propagator. The structure of these 
amplitudes is very simple, as one loop is factorized into a momentum- 
independent expression. 
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The exact analytic computation of the 2-loop electro-weak form factor is a rather 
lengthy process, even though the general method is rather clear. As explained in 
detail in [I], it involves 2 basic steps: 

1. the reduction of the original amplitudes considered before — obtained with 
usual Feynman rules — to a few independent scalar integrals, called master 
integrals (Mis); 

2. the analytic evaluation of all the master integrals generated with the previous 
step with some technique, such as Feynman parameters, dispersion relations, 
small and large momentum expansions or, in our case, differential equations 
in the external kinematical invariants. 

2.1 The reduction to master integrals 

By reduction to master integrals we mean the reduction of all the Feynman diagrams 
to a minimal set of independent scalar integrals. This process involves, in turn, the 
following three steps: 

1. Projection on invariant form factors. 

The amplitude associated to a given Feynman diagram is, in general, a tensor 
integral of the form: 

•^fiuij- / " ft^l" 1^2 p,n2 p,n3 7-)n4 nns ptne ' ^'^ ) 

J -^ii ^ii ^is id 

where /xz/zj . . . denotes collectively spinor indices, four- vector indices and so 
on. For a scalar probe, with an interaction for example of the form ffX, the 
diagram has 2 spinor indices, while for a vector probe, with an interaction 
for instance of the form the diagram has also a 4-vector index. By 

usual form-factor decomposition of Q, one can limit himself to consider scalar 
integrals of the form: 

_ /" tD; jDi N{Si, S2, S3, S4, S5, Sq, Sj) , , 

- a Hid ^2 p,n2 r)"3 r,ni r)"5 n^g l^j 

J il 12 i-i ii is id 

where 

5"! = kl, S2 = kl, S3 = ki-pi, Si = ki- p2, (5) 
85 = k2-pi, Se = k2- P2, Sj = ki- k2. 

The function N{Si, S2, S3, 84^, S^, Sq, S-j) is a polynomial in the kinematical 
invariants S^. 

2. Rotation to independent scalar amplitudes. 

The scalar amplitudes in Eq. (@)) are not linearly independent on each other 
as they involve for instance 6 denominators and 7 scalar products, i.e. a total 
of 13 factors, while there are only 7 kinematical invariants depending on the 
loop momenta ki and ^2- There are two different methods to eliminate the 
redundant amplitudes: 
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(a) Auxiliary diagram or auxiliary denominator scheme jB]. 

We express the scalar products in the numerator of (j3)) in terms of the 
original denominators Di {j = 1 . . . 6) and of an auxiliary denominator 
by means of invertible relations of the form: 

7 

Sk = ^akjDij (6) 
i=i 

Our task is then shifted to evaluate independent scalar amplitudes con- 
taining formally only denominators 

<Sl = I <f'ki(f'k2 ni j^n2 Tjn3 j^m jjns jjn^ jjm 5 (7) 

with rij < 1 for z < 6 and ^7 < 0. Since the numerator is expressed 
by means of ® as a polynomial in the denominators, a term Di in the 
denominator can be canceled by an equal term in the numerator. 
This means that line i is shrunk to a point and the topology number t 
of the diagram is decremented by one®. A given Feynman diagram then 
generates a pyramid of sub-diagrams corresponding to all the possible 
contractions of i > 1 internal lines. All the possible resulting diagrams 
are plotted in Figs. El IH El and El 

Since the contracted propagators may be massless as well as massive, 
an amplitude with for example 2 massive propagators will generally con- 
tain, in its decomposition, independent subamplitudes with 2,1 and 
massive propagators. The amplitudes computed in our previous work 
containing or 1 massive denominators, then, usually appear in the 
decomposition of the present amplitudes; 

(b) Shift scheme [ZIEIE]. 

No auxiliary denominator is introduced. A routing of the amplitudes 
has to be initially assumed with one internal line carrying momentum 
ki and another internal line carrying momentum /c2- This way one can 
simplify powers of k\ and k\ in the numerator against denominators, 
by means of formulas like kf/{kf + a) = 1 — a/{kf + a) with i = 1,2 
and a = 0,m^. A set of cancellation rules between scalar products in 
the numerator and denominators is established. An example of such 
a cancellation is: ki ■ pi/ikf + 2ki ■ pi) = 1/2 — kl/2/{kl + 2ki ■ pi). 
We then make all the possible scalar-product denominator cancellations. 
After such cancellations, the amplitude may loose the denominators with 
momenta ki and ^2, which where initially present by construction. The 
latter are then reproduced with proper shifts of the loop momenta, hence 
the name of the method. In general, with the shifts, new denominators 



^We define the topology number t as the number of different denominators Di occurring in a 
diagram with > 1, irrespective of the powers of the scalar products 
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are generated and the table of scalar-product denominator cancellation 
rules has to be updated to include also the new denominators. We go 
on through steps of cancellations and shifts till no more cancellations 
are feasible. The amplitudes generated in the final step represent the 
independent scalar amplitudes. Since the total number of invariants is 
7, an independent amplitude with t denominators has in general t = 
7 — t scalar products in the numerator. As with the previous scheme, 
it is clear that sub-amplitudes with all the possible contractions of the 
internal lines of the original (linearly dependent) amplitude do occur in 
the decomposition. 

The contractions of internal line(s) in an amplitude may lead to the following 
general simplifications: 

(a) The internal line(s) between 2 external lines are contracted. 

This implies that 2 external lines meet in a new effective vertex and the 
3-point function effectively simplifies to a 2-point function. There are 2 
possibihties: 

• contraction between external fermion lines. 

Since we are computing a 3-point function with 1 general external 
momentum q = Pi + P2 and 2 light-cone momenta pi and p2'- p\ = 
P2 = 0, we obtain a 2-point function with the general momentum q 
flowing through it; 

• contraction between an external fermion line and the probe line. 
The simplification is even greater in this case. The resulting diagram 
depends only on a single light-cone momentum, i.e. on pi or on 
P2-, but not on both. A single light-cone momentum is equivalent 
to a null momentum, as is clearly seen with Wick rotation. The 
resulting diagram is then effectively a vacuum amplitude, in which 
we can set to zero the external light-cone momentum, i.e. pi — or 
P2 — ^ 0. The above property is the basis of many relations between 
amplitudes; 

(b) the internal line where loops overlap is contracted. 

The amplitude factorizes into the product of two 1-loop amplitudes, 
whose computation is trivial. 

3. Reduction of independent amplitudes to master integrals. 

It is possible to reduce the independent amplitudes to a small subset of them 
by means of identities obtained with integrations by parts pUj . 

• In the auxiliary diagram scheme, these identities are obtained as: 







/ 
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with i = 1,2 and v = ki,k2,Pi,P2- By explicitly performing the deriva- 
tives and re-expressing the generated amphtudes in terms of independent 
ones as described before, one obtains relations among independent am- 
plitudes with shifted indices rij ^ rij ±1. 

• In the shift scheme, the identities are obtained as: 

Once the integration- by-parts identities and eventual symmetry relations jHl ^ 
have been generated, the next step is to combine them to achieve the complete 
reductions to the Mis. At present, there are 3 techniques to solve the ibps: 

(a) symbolic method [Til] . 

Historically, this has been the first method used, introduced together with 
the ibps identities themselves. The ibp identities are treated as formal 
recursion equations in the indices nj of the denominators (and of the 
scalar products) to take them to some reference values such as typically 
0,1,2,-1; 

(b) numerical-indices method [Zj. 

Explicit numerical values are replaced for the indices nj. In the auxiliary 
diagram scheme, one takes for the indices of the denominators both posi- 
tive and negative values: rij = ■ ■ ■ — 2, —1, 0, 1, 2 ■ ■ -, while for the indices 
of the auxiliary denominators only negative values: riji = 0,-1, —2 ■ ■ -. 
In the shift scheme, one takes positive values for both the indices of the 
denominators and the scalar products: nj = 0, 1,2---, Ij = 0, 1,2---. 
In either scheme, a homogeneous system of linear equations is then gen- 
erated, whose unknowns are the amplitudes themselves. The system is 
then solved with the method of elimination of variables, after having es- 
tablished the order in which amplitudes and equations have to be solved. 
The amplitudes which remain on the r.h.s., after all equations have been 
used, are the Mis for the specific hierarchy assumed; 

(c) integral method |TT]. 

An integral representation for the expansion coefficients of the indepen- 
dent amplitudes in terms of master integrals is derived. This technique 
has not been used up to now, as far as we know, for original computations 
but for general studies. 

The reduction through the ibps of a given independent amplitude involves Mis 
with subsets of denominators of the starting amplitude. In other words, the 
reduction involves Mis in which 0, 1, 2 ■ ■ ■ internal lines of the original ampli- 
tude have been shrunk to a point. All the Mis appearing in our reductions are 
given in Fig. [7| 
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(a) (b) (c) (d) 




(e) (f) (g) 



Figure 3: The set of 7 independent 6-denominator diagrams, with 2 and 3 massive 
propagators. 

2.2 The evaluation of the master integrals: the differential 
equation method 

In the previous Section we outhned how to reduce the 2-loop Feynman diagrams 
of the electroweak form factor to the master integrals. Let us now sketch how 
to exphcitly evaluate the master integrals with the differential equation method 
[la ESI CH Ca C^. We may identify 2 basic steps: 

1. Generation of the differential equations. 

This step is rather "automatic" in the sense that it does not offer specific 
difficulties. It involves, in turn, the following steps: 

(a) We take the derivatives of the Mis with respect to x at fixed a by dif- 
ferentiating the integrand with respect to the external momenta pi or 
P2; the derivatives of the Mis are then expressed as a superposition of 
linearly-dependent scalar amplitudes; 

(b) We rotate the dependent scalar amplitudes to independent ones, either 
in the auxiliary diagram scheme or in the shift scheme; 

(c) We reduce the independent amplitudes to Mis by using the ibps identities 
and other eventual symmetry relations. This way the system closes on 
the Mis themselves. Linear systems of first-order differential equations 
with variable coefficients are then generated, whose unknowns are the 
Mis themselves. 
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(a) 



(b) 



(c) 



(d) 



(e) 



'"{^^ ^oc> ^co^ "{3^ 



(i) 



(j) 




(k) 



Figure 5: The set of 11 independent 4-denominator diagrams, with 2 and 3 massive 
denominators. 







(a) 



(b) 



(c) 



(d) 



(e) 



Figure 6: The set of 4 independent 3-denominator diagrams together with the single 
2-denominator diagram, with 2 and 3 massive propagators. 
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(a) (b) (c) (d) (e) 




Figure 7: The set of 25 Mis with 2 and 3 massive denominators. Out of them, 4 are 
the product of 1-loop Mis. The dot on a line indicates a square of the corresponding 
denominator. 



2. Solution of the differential equations. 

In order to have a unique solution — as it should — initial conditions can 
often be obtained by studying the behaviour of the Mis close to thresholds 
or pseudothresholds. The system of differential equations is then solved by 
recursion in e. Let us remark that this step is not "automatic", in the sense 
that obtaining a solution is in some cases a matter of intuition and experience. 

Let us make a few comments about specific characteristics of the present computa- 
tion. The most complicated cases involve topologies with more than 1 Ml: 

• t = 4: we encounter 2 vertex topologies having 2 Mis; 

• t = 5: we encounter 2 vertex topologies having 2 Mis and 1 topology having 
3 Mls^ \ 

A convenient basis for topologies with 2 Mis consists of the following amplitudes: 

1. the "basic" amplitude, i.e. the amplitude with unitary numerator and with 
all the denominators having unitary indices; 

2. the amplitude with unitary numerator and with one of the denominators 
squared, or with an independent scalar product left on the numerator. 

The resulting system of 2 differential equations is often triangular in D = A, allowing 
for a simple solution. For the only case of 3 Mis, we found that a convenient basis 
comprises the following amplitudes: 

1. the "basic" amplitude defined above; 

2. the amplitude with unitary numerator and one denominator squared, as in the 
previous case; 

3. the amplitude with an irreducible scalar product in the numerator and the 
denominators all with unitary indices. 

With this choice, the system of 3 ordinary differential equations becomes triangular 
in Z) = 4, allowing for a simple recursive solution in e. 

3 Harmonic Poly logarithms of one variable 

As discussed in the introduction, amplitudes with 2 and 3 massive propagators may 
have thresholds in s = 4m^ and pseudothresholds in s = — 4m^. As a consequence, 
a generalization of the usual harmonic polylogarithms of one variable is needed. 
In the next Section we summarize the salient features of the ordinary harmonic 
polylogarithms (HPLs) that we want to preserve with the generalization, while in 
Section 13.21 we discuss such a generalization. 

^The t = 6 crossed ladder topology also has 3 Mis. 

^We have also used the symbolic method to reduce the topology with 3 Mis, with similar results. 
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3.1 Ordinary Harmonic Poly logarithms 

The basic idea of the HPLs is that of representing a given integral introducing a 
minimal set of transcendental functions by making only linear transformations, such 
as partial fractioning and integration by parts. Furthermore, the transcendental 
functions are defined once and for all. Let us begin with a simple example. The 
integration of the power function with an integer exponent n, 

X x""^^ 1 

x' dx = if n ^ —1 (10) 

n+ 1 n+ 1 ^ ^ 

= logx if n = —1 , 

is again a power function, except for the case n = — 1 for which we have the log- 
arithm. Integration then introduces 1 transcendental function, plus the original 
set of elementary functions. Larger sets of transcendental functions are needed to 
represent double integrals, such as for instance 

dx' r dx' r' dx" 

-logx' = — = -logx log(l -x) -Li2(x), (11) 

1 — x' Jq 1 — X Ji x" 

which involves a new transcendental functions, the well-known dilogarithm: 

U,{x) ^ - / ^ / = y % for |a;| < 1. (12) 

However, many double integrals do not involve new transcendental functions with 
respect to the ones needed to represent single integrals such as (fTUj). Let us indeed 
consider: 

dx' I r dx' r' dx" f I \ 

■logx'= / —I — = log(l-x)+ 1 logx. (13) 



\ — x'Y Jo i^ — x')^Ji x" \1 — X 

The reduction above is done with an integration by parts: we integrate 1/(1 — x)^, 
whose integral is elementary, and differentiate logx, obtaining a simpler, elementary 
function. We then make a partial fractioning of 

111 

(14) 



x(l — x) X (1 — x) ' 

i.e. another linear transformation. In general, whenever we have to integrate the 
product of a transcendental function times an elementary function whose integral 
is also elementary, we do a by-parts integration in order to derive, i.e. to sim- 
plify, the transcendental function. This applies in particular to the integral of the 
transcendental function itself: 

I fix' log x' = X log X — X (15) 
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and in general to all the integrals of the form: 

dx' 



with n ^ 1. The examples given above generalize to multiple repeated integrations, 
of the form: 

dxi dx2 dxi 

with ai ■ ■ ■ a/ some constants. One has to introduce transcendental functions given 
by the repeated integrations of the inverse linear functions + a^): 

dxi dx2 dxi 



Xi + ai J X2 + 0,2 J Xi + ai 

For many applications, such as our previous computation of 2-loop electroweak 
amplitudes with 1 massive propagator [T], it is sufficient to introduce the following 
set of basis functions: 

g{0,x) = -, (19) 

X 

9{l,x) = (20) 
1 — X 

= T^- (21) 

1 + X 

Note that ^'(0, x) has a non-integrable singularity for x ^ 0, while the other functions 
are finite in the same limit. The HPLs of weight 1 are defined as integrals of the 
basis functions: 

r df 

HiO,x) = y - = logx, (22) 
/■^ df 

H{l,x) = j -^ = -log(l-x), (23) 

r dt 

H{-l,x) = -_ = log(l + x), (24) 

Note the slight asjTiimetry related to the lower bound of integration for if(0,x). 

Let us now define the general HPL H{w] x), where w is a vector with w compo- 
nents, consisting of a sequence of "1", "0", and "-1". The weight of a HPL is the 
number of its indices, coinciding with the number of the repeated integrations. 
The HPLs of weight w + 1 have the following integral recursive definition: 

/•X 

H{a,w;x)= f{a;x')H{w,x'), (25) 
Jo 

except for the case with all the indices zero: 

i^(0^;a:) = ^log'"x, (26) 
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or, if a recursive definition is preferred: 

H{0, 0^; x)= [ /(O; x')H{0^, x') = ^ log^+' x, (27) 

where 0^ = (0, 0, • • • , 0) is a vector containing w zeroes. 

3.2 Generalized Harmonic Poly logarithms 

The basis functions for the GHPLs involve various kinds of extensions of the basis 
functions defining the HPLs. 

A trivial extension concerns functions involving different real constants: 

^(-4,^) = (28) 
g{A,x) = (29) 

The above functions are related to amplitudes with threshold/pseudothreshold at 
s = ±4m^ respectively, just as the old functions g{^l', x) represent amplitudes with 
threshold/pseudothrcsold in s = ±m^. We will see in a moment however that the 
above extension in not sufficient. 

To represent 3— point functions with 3 massive propagators one also needs to 
introduce basis functions involving a complex constant:^ 

X 2 * 2 

9ic,x) = _ / ■ (31) 

Note that ^ ± = exp(±i7r/3) are the non trivial cubic roots of and are the 
inverse of each other. 

The non trivial extension, however, involves: 
• radicals of the form 

9{-r,x) = ] , (32) 

y/x{A + X) 

9(1-, x) = —===; (33) 
y/x{4: - x) 

These functions also describe amplitudes with thresholds and pseudothresholds 
in s = ±4m^ respectively. It is indeed well-known that amplitudes involving 
2 particles with the same mass m ^ contain square roots of similar form as 

a phase-space reduction effect; 



^An alternative representation uses only real functions at the price of introducing squares of 
the integration variable: l/(a;^ — a; -|- 1) and x/ {x^ — x + 1). 
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products of radical with inverse linear functions:^° 



g(-l-r,x) = ^ , (34) 

g{l-r,x) = \ -, (35) 

•\/x(4 + x) (1 - x) 

g{-l + r,x) = \ -, (36) 

•\/a;(4 - x) (1 + x) 

+ = r- (37) 

y x(4 — a;) (1 — a;) 

In total, we have added 10 functions to the old basis. All the new basis functions 
above arc finite in x = or have at most an integrable singularity ~ 1/ ^/x for a; 0. 

The GHPLs of weight 1 are given by integrals of the basis functions and can be 
written, in general, in terms of logarithms of complex functions of x: 

r (it 

H{-A,x) = J ^ = log (4 + a;) -2 log 2, (38) 

r dt 

H{A,x) - __^_iog(4-x)+21og2, (39) 



H{c,x) - / 
Jo 

ig (^x - 1/2 - - log (^-1/2 - iV3/2^ , (40) 

H{c,x) = r 

^0 



2 2 

log, 

dt 

t — 2 + 



log (x - 1/2 + iv^/2) - log (-1/2 + iV3/2^ , (41) 



= 21og(V^T4 + ^/^) -21og2, (43) 



rr, ^ r dt ^ . 

H(r, x) = I , = 2 arcsm , (44) 
H(-l -r,x) = / , = arctan W- , 46) 



I 



yi(4 + t)(l + 1) a/3 

V4 + a; — i\/3a; 



Alternative basis functions havi ng the sam e asymptotic hmits ~ 1/a; for x ^ 00 of the previous 
ones can be taken for instance as -\/(4 ± x)/x/{l ± x). 
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H{l-r,x) = / ^ = iog<^^^^^Z^^, (4j 

H{-l + r,x) = I = ^ arctan I W— ^ 1 , (49) 



(50) 



^(1+.^) - I -^__^^log^^^^^ (51) 

Note that all the integrals above have zero as the lower limit of integration. 

The GHPLs of general weight w > 1 are defined exactly in the same way as 
the HPLs, according to Eq. (j^ . where now the components of w can also take the 
values ±4, c, c, ±r, ±1 ± r. 

Let us now discuss how the integrals appearing in the evaluation of the Mis can 
be reduced to the GHPLs plus, of course, elementary functions. As far as the indices 
±4, c and c are concerned, there is basically nothing new with respect to the HPLs: 
the only difference is that individual terms in the expressions for the Mis are in 
general complex and only their sum is, at it should, real. 

When radicals, i.e. semi-integer powers are involved, the situation is more com- 
plicated. The integral of a semi-integer power, 

[ _ I (x + a)""*"^/^ 

/ (x + a)" ^I'^dx = — (n is an integer), (52) 

J n+l/2 

is always a semi-integer power, i.e. an elementary function. There is no need, then, 
to introduce basis functions of this kind. On the other hand, integrals of products 
of radicals, such as 

I{a, p)= l{x + aY{x + hfdx, (53) 



with a and (3 half integers, involve in general new transcendental functions^^. By 
taking a = and h = ±4 (or a = ±4 and 6 = 0), we cover with I{a,(3) the cases 
related to the basis functions with indices ±r. The idea is to shift the indices to 
reference values, such as a = —1/2 and f3 = —1/2 in our (arbitrary) choice of the 
basis functions, by using recursive relations. 

Let us consider the general case in which both indices differ from their target 
values: a 7^ —1/2 and (3 7^ —1/2. The integrand of ()53|) is conveniently rewritten as 

{x + aY{x + hf = ^ {x + aY{x + h)\ (55) 

a/(x + a)(x + h) 



^^Ow the contrary, products of integer powers do not require the introduction of new basis 
functions because partial fractioning completely disentangles the factors: 

1 1 A„ , A. Bn , , Bi 



(a; + a)" (x + 6)'= (x + a)" x + a (a; + &)" 



(54) 
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with n = a + 1/2 and / = /3 + 1/2 general integers. 

We first make an algebraic reduction on the integrand to take one of the two 
indices n and I to its target value zero. There are two possibilities: 

• If 1 of the 2 indices is positive, let's say n > 0, we expand {x + a)" in powers 
oi X + b with the binomial formula: 



s=0 

so that the integrand takes the form 



(56) 



(x + a)°(x + 6)^ = ^ (""Va-fe)"-^ (x + a)-^/2(x + 6)^+'-V2 (jf ^ > q), 

(57) 

in which the index a reached its target value. 

If 6ot/i indices are negative, n < and / < 0, we do ordinary partial fractioning: 
' ' - + + — (58) 



(x + a)l"l (a; + 6)l'l (a; + a)l"l (x + a)l"l-i x + a 

Bill B\i\_i B\ 



SO that 



{x + hf\ {x + hf\-^ x + b' 

{x + ay{x + b)f^ = A|„l(x + a)-Wn|(2; + 6)-V2 
+ >lH_i(x + a)V2-H(^ + 6)-V2 

+ ••• + A, (x + a)-^/^ (x + 6)-V2 
+ S|^|(x + a)-V2(^ + 6)-V2-|^| 

+ S|,|_i(x + a)-V2(x + 6)WI 

+ ••• + B, (x + a)-V2 (2: + 6)-3/2 

(if n< and I < 0). (59) 

In any of the above terms on the l.h.s. one of the two indices reached its target 
value. 

To take also the second index to its target value, we use the integral identity: 

I{a, (3) = -—^{x + ar^\x + bf - ^^^^/(«, - !)• (60) 

The above relation is obtained by doing an integration by parts on I{a,f3): we 
integrate (x + a)" and differentiate (x + 6)^ with respect to x; we then simplify the 
ratio {x + a)/{x + b) as 1 + {a — b)/{x -\- b). 

The first term on the r.h.s. is a finite one and is therefore known. The above equa- 
tion can then be used to lower the index (3 by one unit. For example /(a, 3/2) can 
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be transformed to a linear combination of I{a, 1/2) and I{a, —1/2), while I{a, 1/2) 
can be reduced to /(a, —1/2). In general, by recursively using equation (plj) . one 
can reduce any f3 < —1/2 to (3 = —1/2^^. 

By solving Eq. (jUU]) with respect to the last term on the r.h.s, I{a,P — 1), and 
shifting /? — i> /? + 1, one obtains an identity to raise the index /? by one unit:^^ 

(/? + l)(a-6) (/3+l)(a-fe) '("'^ + ^^- 

By using Eqs. ()60|) and we can then take the index j3 to any desired reference 
value, such as the value —1/2 of our basis. Then, we succeeded in reducing an 
integral of the form to if(±r; x) plus terms containing elementary functions 
only. 

In general, we encounter integrals containing both radicals and GHPLs, of the 
form: 

J{a,(3)= I {x + aY{x + hYH{v,w]x)dx, (62) 



where H{y, w; x) is a GHPL whose first index v has been separated out for a reason 
that will become clear soon. Since w has w components H{v, w; x) has weight w + 1. 
As with the previous integral, we can transform one of the two indices a and j3, let's 
say a, to its target value a = —1/2. The identity to reduce the second index is: 

J(a,/3) = ]- — -{x + a)"+\x + bfH{v,w;x) (63) 

a + p + 1 

1 



a + p + 1 
f3{a - b) 
a + /? + 1 



X + a)"'^^{x + b)^g{v; x)H{w; x)dx 
J(a,/3- 1). 



The second term on the r.h.s. involves the integration of a GHPL of smaller weight. 
We then consider the last term as the only relevant one in the recursion. As before, 
the above equation can be directly used to lower the index (3 by one unit. Note that 
the indices inside the GHPL are not touched by the reduction. Analogously to the 
previous case, by solving Eq. (jU^ with respect to J(a,/5 — 1) and sending (3 ^ (3 + 1, 
one obtains an identity to raise [3 by one unit. 

In some cases, it is also necessary to integrate expressions involving 3 factors, of 
the form: 

L(a, I3]k)= I dx {x + a)" {x + bf {x + c)-^ (64) 



with a and (3 half integers and k a general integer^^. For our computation, the 
relevant cases are a = 0, 6 = ±4 and c = ±1 (or a = ±4, 6 = and c = ±1). 



-^^Since we can take a = —1/2 and the equation is used for /3 > — 1/2, it holds that a + f3 > —1, 
so the singularity m a + (3 — —1 is always avoided. 

^•^Note the the singularity for /? = — 1 is never reached since f3 takes half- integer values only. 

^*We have put a minus sign in front of k just for practical convenience, in order to simplify the 
forthcoming results a little bit. 
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Let us consider the various possibilities for the k index. For k < one can reduce 
the integral ()64j) to the simpler form of I{a,(3) in ()53|) by means of the binomial 
expansion of (x + c)'*^' in powers of a; + a or of a; + 6: 

^'^^ /I jH\ /I jH\ 

(x + c)l^l = Y.[ ) ~ af^-'{x + ay = Y,\ ) (c - by^^~'{x + by. (65) 

s=0 V / s=0 V / 



Therefore we just need to consider the case k > from now on. 

By means of the algebraic relations in (j57|) and (j59|) . we can shift one of the 
half-integer indices a and (3 to the target value —1/2; let us assume for instance 
that a = —1/2. We can then assume the integral of the form: 

L'{P;k) = Li-l/2,P;k) = [ U + bfix + c)-' , (66) 

J ^y{x + a)[x + o) 

with / = /5 + 1/2 a general integer. 

With an algebraic reduction analogous to the one in Eqs. (j57j) and (j59|) . we can 
reduce the product (x + 6)'(x + c)~'^ to a linear combination of terms involving 
powers of either [x + b) or (x + c), but not the product. The terms not containing 
(x + c) belong to the simpler class of the integrals I{a,P) defined in Eq. 
whose reduction has already been discussed. The terms not containing (x + b) have 
the indices a and (3 both equal to their target value —1/2. Our task is then the 
evaluation of integrals of the form: 

m = L(-l/2,-l/2;fc) = / -=i^==(x + c)-^ (67) 

J A/(x + a)(x + o) 

with k integer and strictly positive, k E N . The required identity is: 

L{k)= - {x + a)-^'^{x + b)-^'\x + cy''+^ (68) 



k-l 
1 



2(A;-1) {c-af-^ 
1 1 



+ 



2(A;-1) (c- 6)^-1 
^ fc-i 



.-1 . ^ 1 



2{k 



1=1 



dx (x + a)-3/2 (x + b)-^'^ 
dx (x + a)-i/2 (x + b)-^'^ 

L{1). 



(c — aY~^ (c — by 



The second and the third terms on the r.h.s. do not contain any power of (x + c), 
can be reduced with the previous identities and are therefore known terms. The last 
term contains integrals L{1) of the same form as the one on the l.h.s.: L{k). Eq. (|68p 
is then to be used to relate L integrals, treating the other terms as known functions. 
By using Eq. (jUHj) for instance for = 3 one can reduce L(3) to a superposition of 
1/(2) and 1/(1), while using it for /c = 2 one can reduce 1/(2) to 1/(1). In general. 
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by recursively using Eq. (j68|) . one can reduce any k > 1 to k = 1, i.e. to the basic 
integral^^ 

L(l) = [ dx ^ , (69) 

J y' {x + a){x + b) {x + c) 

which defines the H{±r ± 1; x)'s. 

The derivation of Eq. is analogous to the ones of the previous identities: 
one has to integrate (a; + c)~'^ and to differentiate {x + a)~^^'^(x + b)~^^'^ with respect 
to X. The resulting expressions have to be simplified using {x + c)/{x + a) = 
1 + {c — a)/{x + a) and an analogous equation with b replacing a. One also needs 
the following result: 

1 1 _ ^ 1 1 11 

(x + c)^ x + a~ ^ (c - ay+^ {x + c)''-^ ^ (c - a)'' X + a ^ ^ 

and an analogous equation with b replacing a. 

Finally, the last class of integrals we need to consider is: 

f ^ + + 6/ + (Tl) 

As with the previous case, we can reduce ourselves to the case a = —1/2, /3 = —1/2 
and A; > 0, i.e. to: 

r dr 

F{k) = F(-l/2,-l/2;A;) = / , = H{v,w-x). (72) 

J a/(x + a) (x + 6) (x + c)^ 

The identity reads: 

F{k)= - {x + a)-^/\x + b)-^/\x + c)-''+^H{v,w;x) (73) 



k-l 

1 If 

— — / dx (x + a)-^/2 (x + b)-^^^ H{v,w;x) 

- 1) (c - a)"-^ J 

dx {x + a)"^/^ {x + 6)"^/^ H{v, w; x) 



2{k 

1 



2{k-l) (c-6)'=-i 

+ / dx (x + a)"^/^ (x + 6)"^/^ (x + cV^'^^giv] x)H(w; x) 

k-l J 



^ 2{k 



1 1 



+ 



(c — a)'^ '■ (c — b)'^ 



k-l 



F{1). 



The fourth term on the r.h.s. involves the integration of an elementary function 
times a GHPL of weight w, while the l.h.s involves the integration of a GHPL of 
weight w + 1. This term then has a smaller recursive weight than F and has to be 
considered as a known integral. For the rest, analogous considerations to the one of 
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Note that the singularity of the coefficients in Eq. (|68|l for fc = f forbids any further reduction. 
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the previous identity hold. By recursively using Eq. (f7H|l . one can reduce any k > 1 
to k = 1, i.e. to the basic integral 



which defines the GHPLs of weight w + 2. 

3.3 Closure under the transformation x ^ 1/x — a further 
extension 

The extended basis function set introduced in the previous Section preserves the 
following fundamental properties of the HPLs: 

• the uniqueness of representation as repeated integration noted before 12]; 

• the fulfilling of an algebra. This means that the product of two GHPLs of 
weights wi and W2 can be written as a sum of GHPLs of weight wi + W2. 

The ordinary HPLs, however, have also the property of "closure" under the inverse 
transformation x — >■ 1/x, i.e. anyone of the functions g{±l; x), g{0; x) transforms 
into a linear combination of the same functions under x — > 1/x. This property, 
which is useful for the large momentum expansion of the Mis, is not preserved by 
the generalization discussed in the previous Section. One can impose this property 
on the GHPLs at the price of a second basis extension. Before doing that, however, 
let us discuss in detail how the expansion of the ordinary HPLs for large values of 
the argument x is performed. The foUowings steps are taken in order: 

1. we set x = l/y in the HPL under consideration: 



2. we use the identities which allow to reduce H{w; 1/y) to a combination of 



3. we expand H{w'; y) for y — > 0, i.e. for small value of the argument, as explained 
in some detail in Section IHl 

4. we perform the inverse substitution y = 1/x in the final result. 
As an example consider the simple weight- 1 HPL 




(74) 



H{w;x) = H{w; 1/y) ; 



(75) 



H{w';yys- 




(76) 



We have: 



H{-l;x) 



i-A/y) 



(77) 




(79) 



(78) 



HiO-y)+Hi-l-y), 



(80) 
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where, moving from Eq. (j78|) to Eq. (j79p . we divided the integral into the sum of 
two integrals (the first between and 1 and the second between 1 and 1/?/) and we 
replaced t by 1/t'. 

Expanding the r.h.s. of Eq. ()80|) in series of y and re-expressing y as 1/x we 
finally have: 

H{-l;x) "=~logx+---^ + of^ 



f82) 



X 2x^ \ x^ 

For the weight-2 HPL 



HiO,-l;x) = J^^jHi-l;t) 



we have: 



H{0,-l-x) = HiO,-l;l/y) (83) 

/"^/^ dt 

= ff(0, -!;!) + J jH{-l;t), (84) 

n dt' 

= if(0,-l;l)- / —Hi-l;l/t'), (85) 



^6) 



= //(0,-l;l)- l'^[-H{0;t')+Hi-l-t')], 
= 2H{0,-l;l) + H{0,0;y)-H{0,-l;y), (87) 

where, moving from Eq. (j85|) to Eq. (j86|) we used Eq. (j80|) . 

Expanding the r.h.s. of Eq. ()87|1 in series of y and re-expressing ?/ as 1/x we have 
finally: 



H{0, -1; X) = 2H{0, -1; 1) + - log^ x-- + — + 0^- 

Let us remark that performing the substitution x = l/y the basis functions 
g{n, l/y) are expressed in terms of the basis functions belonging to the same set, 
g{n',y): this is the closure under the transformation x 1/x. Moreover, in the 
final result the constants H{w; 1) do appear. The latter have been systematically 
evaluated and tabulated in jTTj up to weight 4 included. 

Following exactly the same steps in the case of the GHPLs, we find that the set 
of basis functions given in the previous Section is too small to preserve the closure 
under the transformation x — > 1/x. We can understand it, for example, trying to 
expand the GHPL 

r dt 

HiA; x) = __ = - log (4 - x) - 2 log2 . (89) 

We have: 

H{A-l/y) = ff(4; 1) + ^ , (90) 
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y dt' 1 



^(4; i; 

/y 
dt' 



4 t' t'-\' 



V t' 



H{A-l) + H{Q-y) 



y dt 
dt' ^ 



4 



(91) 
(92) 
(93) 



The integral in Eq. does not belong to the set of GHPLs of the previous Section. 
We then add to the basis functions defined in Eqs. (j2HH2ZI) the following functions: 



g{±l/A-x) 
g{±l ± r/4; x) 

fif(ro/4;x) 

5((-ro/4;x) 



1 



=F X 



1 - 2^^ X 



X \ I X ^ 



l-2jx+\ 



(94) 
(95) 

(96) 

(97) 



X\ X + 



To avoid a pole in x = 0, we have subtracted (^(0; x) = 1/x with a proper coefficient, 
the residue in a; = of the new functions. In taking the limit x — > inside square 
roots one has to remember that x x — ie, with e a small positive number. The 
related GHPLs of weight 1 are defined as usual: 



H(w: x) 



g{w]t)dt, 



(9J 



since all the functions (jSUHZI) are integrable in x = 0. 
As an example, the reader may verify that: 

H{4,-r;x) = H{4,-r;l/y) 

= -i7(4;l)if(r/4;l) + ^iJ(-l,-ro/4;l) 



(99) 



+i7(l,r/4;l)--/7(-ro/4,4;l) 



+ 



H{r/A; 1) + -H{-ro/4; 1) 



H{0;y) + HiA;y) 



--H{4, -ro/4; y) - -if(-l, -ro/4; y) . 



(100) 



Expanding the GHPLs for ?/ — > and substituting in the final expression y = 1/x 
we finally obtain: 



if (4, -r;x) 



-H{4- l)H{r/4- 1) + Ih{-1, -ro/4; 1) 
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+i7(l,r/4;l)--i7(-ro/4,4;l) 



H{r/A- 1) + -H{-r,/A- 1) 



logx + 0[- 

X 



(101) 



As in the case of the HPLs, the final expression contains the constants H{w] 1), 
which have to be evaluated. The x = 1 GHPLs of weight one, H{a; 1), can be 
expressed in terms of known transcendental constants in an elementary way. In 
Eq. (llOlj) . the coefficient of the leading logarithm is, for example: 



-ro/4; 1) 



TT 



3 -2[log(l 



VS) -log2]. 



(102) 



Similar reductions are also possible in other cases but, in general, it is non trivial 
to reduce the H{w] l)'s with a general weight w to a. minimal set containing known 
transcendental constants such as ({n) and eventual new transcendental constants. 

The reduction of the H{w] l)'s to a minimal set of transcendental constants is 
beyond the goal of the present paper [TS]; we restrict ourselves to the evaluation of 
the asymptotic expansions of the diagrams involving only ordinary HPLs. 



4 Results for the master integrals 

In this Section we present the results of our computation of the Mis involving up to 
6 denominators included, which constitute a necessary input for the calculation of 
the 2-loop vertex diagrams in Figs. ^ and El They are expanded in a Laurent series 
in 

e = 2-L)/2, (103) 

up to the required order in e. The coefficients of the series are expressed in terms 
of GHPL's (see Section El) of the variable x, defined as : 

x = —, (104) 
a 

where s = —q^ is the cm. energy squared^^. It holds q = pi + P2 and we defined 
a = m?. We denote by /i the mass scale of the Dimensional Regularization (DR) — 
the so-called unit of mass. We work in Minkowski space and we normalize the loop 
measures as: 

= . .j^? = . (105) 



fr(3-f) ^7r2-T(l + e) 



This definition makes the expression of the 1-loop tadpole — the simplest of all 
loop diagrams — particularly simple pp. The denominators appearing in the master 
integrals are listed below: 
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We define the scalar product of two 4-vectors as: a ■ b = —aobo + a ■ b. 
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A; 1 ■ 


(106) 








(107) 


T>3 




{ki + k2f , 


(108) 


T>4 





(Pi - ^1)' , 


(109) 


T>5 





(P2 + ^1)^ , 


(110) 






(P2 - k2f , 


(111) 


V7 




(pi — ki + /co)^ , 


(112) 






(p2 + ki — k2Y , 


(113) 






(Pl + P2 — /^l)^ , 


(114) 







fPl + P2 — ^2)^ 7 


(115) 






(Pl + P2 — ^1 — ^2)^ , 

Vi -L X £j f J 


(116) 







/c? + a , 


(117) 






fco + a , 


(118) 






+ k2f + a, 


(119) 








(120) 






(P2 + + a , 


(121) 


Vn 




(P2 - k2)^ + a , 


(122) 






(Pi - /ci + /c2)^ + a , 


(123) 






(P2 + fci - ^2)^ + a , 


(124) 


T>20 




(pi +P2 - /ci)^ + a, 


(125) 


V21 




(Pi +P2 - k2)^ + a, 


(126) 


V22 




(pi +P2 - A;i - k2)^ + a. 


(127) 


rding 


■ : the increasing number of the denominators, which 



corresponds more or less to the level of complexity. 
4.1 Topology i = 3 



= (v)'"E^-';'"+o(e=), 



13^14 



(128) 
(129) 



where: 



7(1) 
-2 



3 

2 ' 



(130) 
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.(1) 



a 
(1) 



— - V^H{r, 0; 1) 



.(2) 



9 
2 

21 

~2 

45 

"Y 

93 
"Y 



V3 



3H{r, 0; 1) + H{r, 0, 0; 1) + if (4, r, 0; 1) 



7H{r, 0; 1) + 3/f(r, 0, 0; 1) + 3/7(4, r, 0; 1) 



ii(r, 0, 0, 0; 1) + H{A, r, 0, 0; 1) + ff(4, 4, r, 0; 1) 



(131) 
(132) 
(133) 

(134) 



This diagram has been originally evaluated in ref. [Tni . 

Even though the above MI is a vacuum amplitude, its computation is non trivial 
because of the presence of 3 massive propagators. We have evaluated it with the 
following method. We consider a vacuum sunrise with 2 internal lines with equal 
mass m and the third internal line with the different mass (squared) m'^ = zm'^. 
We then differentiate the vacuum sunrise with respect to z and rewrite the result 
in terms of Mis by using the ibps identities. The resulting differential equation 
represents the evolution in one of the masses and is solved as in usual cases. We set 
at the end z = 1 to obtain the equal-mass case. 

Let us make a few remarks about the above result: 

• The finite part of the MI, i.e. the 0{e^), involves 1 transcendental constant: 
H{r,0;l), related to the Clausen function; 

• the 0(e) part involves 2 new independent transcendental constants: H{r, 0, 0; 1) 
and H{4:,r,0; 1); 

• the O(e^) part involves 3 new transcendental constants: H{r, 0, 0, 0; 1), if (4, r, 0, 0; 1) 
and ii(4,4,r,0;l). 

To simplify the above expressions, we have used the following identities to move the 
"0" indices to the right: 



H{0,r,0;l) 
H{0,r,0,0;l) 
iJ(0,0,r,0;l) 
if(0,4,r,0;l) 



-i7(r,0;l), 
-2ii(r,0,0;l), 
-3/i(r,0,0,0;l), 
3if(r,0,0,0;l), 

-//(4,0,r,0;l)-2i7(4,r,0,0;l). 



(135) 



The above relations are obtained by transforming products of if 's into linear com- 
binations of if's, as for example in: 



= H{0; l)H{r; 1) = H{0, r; 1) + H{r, 0; 1). 



(136) 
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Topology t = 4 



2{4-D) 



1 



^ i=-2 



(137) 
(138) 



1 

2' 

5 X + 4 



Hi-r;x), 



2 ^/x{x~+4) 

IQ 1 r + 4 



(139) 
(140) 



4i7(_r;x)--i/(-r,-l;x) 



-H{-A,-r;x) 





"1 




- + 1 







i/(-l;a;), 



(141) 



65 1 

- + 5C(2) - C(3) - 3H{0, -1; x) + 2i/(0, -1, -1; x) - -H{Q, 0, -1; x) 



+ 



X 



+ 1 



\^-7H{-l;x)+4H{-l,-l;x)-H{0,-l;x)^ 



x + A 



^/x{x + 4) 



-2(6+C(2))//(-r; x) +6/i'(-r, -1; x)-6H{-r-l-l;x) 



+3H{-r, 0, -1; x) + AH{-A, -r; x) - ^H{-A, -r, -1; x) 



-i7(-4,-4,-r;a;)j>, 

211 Q r 1 

— + 19C(2) + ^C'(2) - 5C(3) - [l3 + C(2)] H{0, -1; x) 

-SH{0, -1, -1, -1; x) + 3H{0, -1, 0, -1; x) + 12H{0, -1, -l;a;) 
-31^(0, 0, -1; x) + 2i/(0, 0, -1, -1; x) - ^/^(O, 0, 0, -1; x) 

I (2(C(3)-4C(2)-16)ff(-r;x) + (6 + C(2))[3iJ(-r,-l;x) 

^x(a; + 4) I 

+2i7(-4, -r; a;)] - 24/7(-r, -1, -1; x) + 24:H{~r, -1, -1, -1; x) 
-9H{-r, -1, 0, -1; x) + 12H{-r, 0, -1; x) - 12H{-r, 0, -1, -1; x) 
+3i/(-r, 0, 0, -1; x) - 6H{-4:, -r, -1; x) + 6H{-4:, -r, -1, -1; x) 

-3//(-4, -r, 0, -1; x) - 4/7(-4, -4, -r; x) + ^//(-4, -4, -r, -1; x) 



(142) 



+ii"(-4, -4, -4, -r;a;) ^ - 



)- 






_x 



{(33 + 2C(2))i/(-l;x) 
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-28H{-1, -1; x) + 16H{-1, -l,~l;x)- 6H{-1, 0, -1; x) 

+7H(0, -l;x)- 4H{0, -l,-l;x) + H{0, 0, -1; x) } . (143) 

The above 2-point function is the simplest amphtude having thresholds both in 
s — and s — Arm?. It does not have any pseudo-threshold. The indices ap- 
pearing in the GHPLs are indeed only 0, —1, —4 and — r. Note that the index — r 
appears eventually only once in the iJ's. The related terms have coefficients always 
containing radicals, in order to reproduce the right causality structure. 



where: 



Z7(3) 



1 

2' 
3 

2' 

^ + //(-I, x) - H{0, -1, x) + - H{-1, x) - C(2)i/(1, x) 



X 



+2H{l,0,-l,x) 



(144) 
(145) 

(146) 
(147) 

(148) 



-- - C(2)+7//(-l, x) - C(2)i/(1, x) - 4//(-l, -1, x) - 2//(0, -1, x) 

+4:H{0,-l,-l,x) + H{0, 0,-1, x) + 2H{1, 0,-1, x) + -\7H{-l,x) 
-(C(2) - C(3))i^(l, x) - AH{-1, -1, x)+H{0, -1, x) - a2)H{0, 1, x) 
+2H{0, 1, 0, -1, x) + 2H{1, 0, -1, x) - 8H{1, 0, -1, -1, x)] . (149) 



where: 




-16^' 



(150) 



(152) 
(153) 
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.(4) 



^1^ 



(4) 



1 + ^H(-1,x)-\h(0,-1,x) 



, (154) 



^ + ^Ci2)-3H{-l,x)+2H{-l,-l,x)-H{0,-l,x)-^ 



^i^(-M) 



+^C(2)//(1, x) - Ih{-1, -1, x) + ^H{0, -1, x) - ^H{1, 0, -1, x) 



+x 
1 



lOQ Q 41 1 Q 

"61 + 8^^^^ " 16^^"^' + 4^(2)i/(l, x) + -H{-1, -1, x) 



.^H{0-l,x)-H{0-l-l,x)-^H{0,0-l,x)-^H{l,0-l,x) , (155) 

SQ 1 25 1 

- + C(2) - -C(3) - -H{-1, X) - Cmi-1, x) - -a2)H{l, x) 

+12//(-l,-l, x) -^/7(0,-l, x) -8i/(-l -1,-1, x)+3//(-l, 0,-1, x) 



+4/^(0, -1, -1, x) - ^ii"(0, 0, -1, x) + /^(l, 0, -1, x) - ^ 



21 
32 



i/(-l,x) 



+^c(2)i/(-i, x) + ic(2)i/(i, x) - ^) - \h{-i, -1, x) 

^H{0, -1, x)+^C(2)//(0, 1, x)+2H{-l-l-l, x) - ^H{-1, 0,-1, x) 
-^H{0-l-l,x) + ^H{0,0-l,x)-Hil,0-l,x)-^H{0,l,0-l,x) 



+2H{l,0,-l,-l,x) 



+ X 



1681 41 , , 3,, , 379 ,,, 
^ + ^C(2)--C(3)--i/(-l,.) 



-^C(2)i/(-l, x) + C(2)//(l, x) - ^C(3)i/(1, x) + ^H(-l, -1, x) 
-^i/(0,-l,a:) + ^C(2)i/(0,-l,a;)-ic(2)i/(0,l,a;)-6if(-l,-l,-l,a;) 
^i/(-l, 0,-1, x)-H{0-l-l, x) - ^-H{0, 0,-1, x)-2H{l, 0, -1, x) 
+4i/(0, -1, -1, -1, a;) - Ih{0, -1, 0, -1, x) + H{0, 0, -1, -1, a;) 
+ ^i/(0,0,0,-l,a;) + Ji/(0, 1,0,-1, a;) + 2i/(l, 0,-1,-1, x) . (156) 



The above 2 Mis, which contain 2 massive denominators, can be expressed in terms 
of ordinary HPLs. The reason is that the 2 massive hnes, roughly speaking, are 
in different channels: one is in the s channel while the other is in the t channel. 
The amplitudes do not have thresholds/pseudothresholds in s = ±4m^, but only in 
s — ±m^. Both the indices "1" and "—1" do indeed appear inside the HPLs. The 
presence of 2 massive denominators is then a necessary but not a sufficient condition 
in order to have thresholds or pseudothresholds in s = ±4m^. 



33 




V ' i— — ') 



where: 



-in — 



^ :^-5C(2)-2C(3)-19i/(0,x) + C(2)i/(0,x)-i/(r,r,0,x) 



1 

2 ' 

19 4 — a; 

- - C(2) - 5H{0, x) + H{0, 0, a;) + —===H{r, 0, x) 
2 ^x(4 - x) 

+-H{r,r,0,x), 

X 

65 



(157) 
(158) 

(159) 
(160) 

(161) 



+5H{0,0,x) - H{0,0,0,x) 



4- X 



V2^(4 - x) 
-H{r, 0, 0, x) - H{0, r, 0, x) + 2/7(4, r, 0, x 



C(2)i/(r,a;) + 5if(r,0,a;) 



X 



a2)H{r,r,x) 



+3H{r, r, 0, x) - H{r, r, 0, 0, x) - H{r, 0, r, 0, x) + 2i7(r, 4, r, 0, x) 
+i/(0,r,r,0,a;)j . (162) 

The double and the simple poles in the MI above have ultraviolet origin. The 
amplitude has indeed an ultraviolet sub- divergence related to the integration of the 
bubble together with an over-all UV divergence. 



(163) 
(164) 




= (^)"e^*'^''' + o(^'). 



where: 








_ 


1 






a 


X 


1 




a 


-H(0, 

X 


x) , 


Fi" _ 
a 


1 

X 


4 


-C(2) 



A- X 



X^/x{4: — x) 



(165) 
(166) 

H{r,0,x), (167) 
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(6) 



2C(2) + 2C(3) + 4H{0, x) - Ci2)H{0, x) - 2H{0, 0, x) + H{r, r, 0, x) 
A- X 



+H{0,0,0,x) ^ 

X^yx[A — X) 

-H{0,r,0,x) + 2H{4:,r,0,x) 



C{2)H{r, x) + 2H{r, 0, x) - H{r, 0, 0, x) 



(168) 



X 



16 - 4C(2) - ^^(2) - 4C(3) - 2(4 - C(2) - C(3))i^(0, x) 



+ (4 - ({2))H{0, 0, x) - Ci2)H{r, r, x) - 2H{r, r, 0, x) + if (r, r, 0, 0, x) 
+H{r, 0, r, 0, x) - 2H{r, 4, r, 0, x) - H{0, r, r, 0, x) - 2H{0, 0, 0, x) 

+H{0, 0, 0, 0, x)] f' -" \2C{2)H{r, x) + 2C{3)H{r, x) +4H{r, 0, x) 

-C(2)i/(r, 0, x) - C(2)if(0, r, x) + 2C(2)i7(4, r, x) + 3i/(r, r, r, 0, x) 
-2i/(r, 0, 0, x) + 4i7(4, r, 0, x) - 2if (0, r, 0, x) + i/(r, 0, 0, 0, x) 
+H{0, r, 0, 0, x) + H{0, 0, r, 0, x) - 2/7(0, 4, r, 0, x) - 2H{4, r, 0, 0, x) 



-2H{4:, 0, r, 0, x) + 4i7(4, 4, r, 0, x) 



(169) 



The double pole in e in the MI above is the product of a simple UV pole coming 
from the nested bubble and of an IR pole coming from the massless line squared. 

The above 2 Mis have a pseudothreshold in s = —Am? related to the exchange 
of 2 massive particles in the t channel and consequently only GHPLs with indices 
0, 4, and r do appear in the e expansion. These Mis have also been computed in 
by means of a transformation well-known in QED eliminating the square roots: 



'l + z) 



X 



(170) 



In general, this change of variable is very convenient for amplitudes not having the 
pseudothreshold in s = — m^. 

4.3 Topology t = 5 



1 



V2V3VUVUV20 



^ ^ i=0 



where: 



aF, 



(7) 



H{0, 0, -1, x) + 3H{-r, -r, -1, x) - 2H{-r, -r, 0, x) 



-2H{-r,0,-r,x) 



(171) 
(172) 



173) 
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aFl 



(7) 



- \2H{0, 0, -1, x) - 6C{2)H{-r, -r, x) + 6H{-r, -r, -1, x) 
X i 

-AH{-r, -r, 0, x) - AH{-r, 0, -r, x) - 4/7(0, 0, -1, -1, a;) 
+H{0, 0, 0, -1, x) - 12H{-r, -r, -1, -1, x) + 6H{-r, -r, 0, -1, x) 
+2H{-r, -r, 0, 0, x) - 3H{-r, -4, -r, -1, x) + 2H{-r, -4, -r, 0, x) 
+2H{-r, -4, 0, -r, x) + 2H{-r, 0, -r, 0, x) + 2H{-r, 0, -4, -r, x) 



+2i7(-r,0,0, -r,x) 



(174) 



The above 2-point function has thresholds in s = 0,m^ as well as in s = 4m^ and 
no pseudothresholds. The index — r always appears twice in the if's and the related 
coefficients contain no radicals. The finite part 0(e°) of F^''^ has been computed 
by the authors of by means of a resummed small momentum expansion. With 
the help of the Mathematica [221 we have computed the first 20 terms of the small 
momentum expansion of our result and compared with their result, finding complete 
agreement. 



1 



2\ 2€ 



(8) 



(175) 
(176) 



where: 



(8) 



C(2) (^H{0, c, x) + 2H{0, 1, x) + H{0, c, x)^ - 3/7(0, c, 0, -1, x) 



-3H{0, c, 0, -1, x) - 4:H{0, 1, 0, -1, x) + -H{r, 0, 1) H{0, c, x) 



X 



-i7(0,c,x) 



(177) 



where: 



{P2 ■ h) 



P2 ■ h 



(178) 



(179) 



.(9) 



-1 + H{-1, x) - C(2)i/(1, x) - -C(2) \H{c, x) + i7(c, x)\ 

-H{0, -1, x) + 2if(l, 0, -1, x) + ^[H(c, 0, -1, x) + if(c, 0, -1, x)] 



+ U\[V3H{r,0,l) + a2) 



H{c,x) + H(c,x) 



-H{c,0,-l,x) 
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--H{c, 0, -1, x) + H{-1, x) + C(2)//(l, x) - 2H{1, 0, -1, x) 



+UHir,0,l)(2-^) H{c,x)-H{c,x) 



+ 



X 



C{2) (^H{c,x) 



-H{c,x)^ -3H{c,Q,-l,x) + 3H{c,Q,-l,x) |, (180) 

-7 - C(2) + 6H{-1, x) - 2C(2)i/(l, x) + C(3)i/(1, x) 
~ (C(2) - C(3)) [h{c, x) + Hie, x)] - 4H{-1, -1, x) 

-HiO, -l,x)- 3C(2)//(0, l,x)- ^C(2) [niO, c, x) + //(O, c, x)] 

+C(2) \h{c, c, x) + H{c,c, x)\ + ic(2) [h{c, c, x) + i/(c,c, x)] 

+2C(2) [//(c, 1, x) + H(c, 1, x)] + AH(0, -1, -1, x) + H(0, 0, -1, x) 



+4i/(l, 0, -1, a;) + - i/(c, 0, -1, x) + i/(c, 0, -1, 



-3 



i?(c, c, 0, -1, x) + if(c, c, 0, -1, x) 



3 

2 L 



-f/"(c, c, 0, —1, x) 



+i7(c, c, 0, —1, x 
1 



-6 



//(c, 0, -1, -1, x) + E{c, 0, -1, -1, a;) 



i/(c, 0, 0, -1, x) + //(c, 0, 0, -1, X 

9 



-4 



i/(c, l,0,-l,x) 



+i/(c,l,0,-l,x) 



+ 



i/(0, c, 0, -1, x) + i/(0, c, 0, -1, 
1 



+6//(0, 1, 0, -1, x) - 8i/(l, 0, -1, -1, x) + [^(r, 0, 1) 

-ii"(r, 0, 0, 1)1 \e[c, x) + (c, x)] + ii"(r, 0, 1) \e[c, c, x) 
+i?(c, c, x) - 2if(c, c, x) - 2if(c, c, x) + if (0, c, x) + i?(0, c, x) 
-i/(0, r, 0, 1) i/(c, x) + E{-c, x) + i/(4, r, 0, 1) i/(c, x) + Ei^c, x)] } 



+ 2 



ii(r, 0, 1) + E{r, 0, 0, 1) + i/(4, r, 0, 1)] [//(c, x) - E{c, x)j 



+3E{0, c,x) - 3E{0,c,x) 



E{r, 0,1) + E{r,0,0,l) 



+E{r, 0, 1) E{c, c, x) - E{c, c, x) - 2E{c, c, x) + 2E{c, c, x) 

1 

2x 

+E{4:, r, 0, 1) E{c, x) - E{c, x)] + E{r, 0, 1) \e{c, c, x) 
-E{c, c, x) - 2E{c, c, x) + 2E{c, c, x) + E{0, c, x) - E{0, c, x)] 

2x 

-C{2)E{0, c, x) - 4C(2)//(c, l,x) + 4C(2)//(c, 1, x) - 2C{2)E{c, c, x) 



(C(2) - C(3)) (^E{c, x) - E{c, x)) + amiO, c, x) 
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+2C(2)iJ (c, c, x) - C(2)//(c, c, x) + C(2)if (c, c, x) - 3H{c, 0, -1, x) 
+3H{c, 0, -1, x) + 8//(c, 1, 0, -l,x)- 8H{c, 1, 0, -l,x) 
-3H{c, c, 0, -l,x) + 3H{c, c, 0, -l,x) + 3H{c, c, 0,-1, x) 
-3H{c,c,0,-l,x) + 6H{c,0,-l,-l,x) - 6H{c,0, -1, ~l,x) 
+H{c, 0, 0, -1, x) - H{c, 0, 0, -1, x) - 3H{0, c, 0, -1, x) 

+3H{0,c,0,-l,x)\\ . (181) 



where: 



14 



(182) 
(183) 



lS^ylH{r,0,l) H{c,x) + H{c,x) 



-3H{c, 0, -1, x) + 3i/(c, 0, -1, x) 



+ i- 



.V3 



'a2)(^H{c,x)-H{c,x) 
(184) 



1 rV3 

'x\ 3 



if(r, 0, 0, 1) + H{4:, r, 0, 1) H{c, x) + H{c, x) 



+H{r, 0, 1) H{c, c, x) + H{c, c, x) - 2H{c, c, x) - 2H{c, c, x) 

+H{0, c, x) + H{0, c, x)] - [C(3) (^H{c, x) - H{c, x)) 

-C(2) (h(0, c, x) - H(0, c, x) - AH(c, 1, x) + AH(c, 1, x) - 2H(c, c, x) 

+2H{c, c, x) - H{c, c, x) + H{c, c, x)) - 2H{c, c, 0, -1, x) 



+2H{c, c, 0, -1, x) - H{c, c, 0, -1, x) + H{c, c, 0, -1, x) 
-AH{c, 0, -1, -1, x) + AH{c, 0, -1, -1, x) - \h{c, 0, 0, -1, x) 

+ Ih{c, 0, 0, -1, x) - ^i/(c, 1, 0, -1, x) + ^H{c, 1, 0, -1, x) 
+H{0, c, 0, -1, x) - H{0, c, 0, -1, x) 



(185) 



The above 3 Mis are real as complex H{- • • c ■ ■ • ; a;)'s always appear in the combi- 
nations: H{- ■ -c - ■ ■■,x) + iJ(- ■ ■ c ■ ■ ■ ; x) and i[H{- ■ ■ c ■ ■ ■ ; x) — H{- ■ ■ c ■ ■ ■ ; x)]. As 
discussed in the previous Section, the above topology is the only one having 3 Mis. 
The amplitudes have a threshold in s = in agreement with Cutkowsky rule as 
well as a pseudo thresholds in s = —m^. 
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where: 



= — //(l,0,x), 



(186) 
(187) 

(188) 



aF, 



(11) 



(11) 



--|C(2)i/(l, x) + 2/^(1, 0, x) - H{r, r, 0, x) + H{0, 1, 0, x) 
-//(I, 0, 0, x) + H{1, 1, 0, x) - 3H{1 + r, r, 0, x) 

+V3H{r,0,l)H{l,x)Y (189) 

i|-2C(2)i/(l, x) - 2C(3)i/(l, ,t) - C(2)if (0, 1, x) - 4H{1, 0, x) 
X L 

+C(2)if(l, 0, x) - C(2)i^(l, 1, x) + amir, r, x) + 3C(2)i/(l + r, r, x) 
-2i/(0, 1, 0, x) + 2/7(1, 0, 0, x) - 2H{1, 1, 0, x) + 2iJ(r, r, 0, x) 
+6i/(l + r, r, 0, x) - i/(0, 0, 1, 0, x) + i/(0, L 0, 0, x) ~ i/(0, 1,1,0, x) 
-i/(l, 0, 0, 0, x) - i/(l, 0, 1, 0, x) + i/(l, 1, 0, 0, x) - i/(l, 1, 1, 0, x) 
-if(r, r, 0, 0, x) - H{r, 0, r, 0, x) + 2iJ(r, 4, r, 0, x) + i7(0, r, r, 0, x) 
-?,H{1 + r, r, 0, 0, x) - 3/J(l + r, 0, r, 0, x) + QH{1 + r, 4, r, 0, x) 

+3i^(0, 1 + r, r, 0, x) + 3/^(1, 1 + r, r, 0, x) - ^3 (^2ii"(r, 0, 1) 

+//(r, 0, 0, 1) + i/(4, r, 0, 1))//(1, x) + H{r, 0, 1) (i/(0, 1, x) 

+i^(l,l,x))]}. (190) 

The above amphtude has pscTidothrcsholds in s = —rn? and s = — 4m^ and is the 
only one containing GHPLs with the index 1 + r. The latter is then related to the 
(virtual) transition of a particle with mass m ^ into a pair of particles with the 
same mass, i.e. a bubble with 2 equal mass lines. The index r appears in the GHPLs 
only or 2 times. 




- (^^fFr + Oie), 



where: 



aFr 



(12) 



X 



C(2)//(0,l,x) + //(0,l,l,0,x) 
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(191) 
(192) 

(193) 




(P2 • h) 



2(4-D) 



P2 ■ h 



(194) 
(195) 



where: 



1 

+ - 

X 



1 + C(2) + H{0, x) + C(2)//(l, x) + H{1, 0, x) + H{1, 1, 0, x) 
a2)H{l,x) - H{l,0,x) - H{l,l,0,x) 



(196) 



(13) 



-7 + 4C(2) - C(3) + 6H{0, x) + 3C(2)i/(l, x) - C(3)i/(1, x) 

-2H{0, 0, x) + 3C(2)iJ(0, 1, x) + 3iJ(l, 0, x) + 5C(2)i/(l, 1, x) 
+iJ(0, l,0,x)- 2H{1, 0, 0, x) + 3iJ(l, 1, 0, x) + H{1, 0, 1, 0, x) 

-2H{1, 1, 0, 0, x) + 5iJ(l, 1, 1, 0, x) + 3iy(0, 1, 1, 0,x) + - \c{3)H(l, x) 

X L 

-3C(2)i7(l, x) - ({2)H{0, l,x)- 3H{1, 0, x) - 5C(2)/7(1, 1, x) 
-//(0, 1, 0, x) + 2H{1, 0, 0, x) - 3H{1, 1, 0, x) - H{0, 1, 1, 0, x) 

-H{1, 0, 1, 0, x) + 2//(l, 1, 0, 0, x) - 5H{1, 1, 1, 0, x)] . (197) 

The above 2 Mis contain HPLs with indices "0" and "1" only. They represent the 
emission of a photon by a charged vector boson in the t channel, and therefore have 
only a pseudothreshold in s = — m^. They arc IR (as well as UV) finite because the 
photon is emitted internally to the basic 1-loop triangle. 




J ViVsVe 



a J 



(198) 
(199) 



where: 



aF, 



(14) 



C(2)//(0, -1, x) + C(2)//(0, 1, x) + H{0, -1, 0, -1, x) 



-2H{Q,l,0,-l,x) 




(pi ■ k2) 



Pi ■ h 



^ ^ 1=0 



(200) 

(201) 
(202) 
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3 + C(2) H{-1, x) + ^C(2)//(l, x) - ^i/(0, -1, x) 



3 1 
-2 + 2 

+\h{-1, 0, -1, x) - H{1, 0, -1, ^) + ^ I 2 



1 1 1 



3 + C(2) H{-l,x) 



-^amihx) + Ih{-1,0, -1,x) + H{1,0, -l,x) } , 



(203) 



-9 - C(2) + I 



15-C(3) 



H{-l,x) 
3 



C(2) + ic(3) 



6 + C(2)] H{-1, -l,x) + -C(2) [i/(0, -l,x) + H{0, 1, x)] 
+6/^(0, -1, -1, x) + ^//(O, 0, -1, x) + 2H{1, 0, -1, x) 
-H{-1, -1, 0, -1, x) - 2i/(-l, 0, -1, -1, - Im-l, 0, 0, -1, x) 
+^H{0, -1, 0, -1, x) - 3H{0, 1, 0, -1, x) + 4//(l, 0, -1, -1, x) 



+ 



1 r 1 



X 2 



15-C(3) E{-\,x)^ 



C(2) + -C(3) 



6 + C(2)] E{-\ -1, 2;) + ^ [3 + C(2)] i/(0, -1, x) - ^C(2)i/(0, 1, 
-2i/(l, 0, -1, x) - -1, 0, -1, x) - 2iJ(-l, 0, -1, -1, x) 

0, 0, -1, x) + ^i/(0, -1, 0, -1, x) + i/(0, 1, 0, -1, x) 



-4/7(1, o,-i,-i,x; 



(204) 




^ ' i=—l 



(205) 
(206) 



-H{—r, —r, x) , 
4//(-r, -r, x) - 3//(-r, -r, -1, x) - 2i/(-r, -4, -r, x) 



+2if(0, -r, -r,x) 
1 



(207) 



(208) 



4(^2 + C(2))i/(-r, -r, x) - 6//(-r, -r, -1, x) + 4//(0, -r, -r, x) 
4//(-r, -4, -r, x) + 12H{-r, -r, -1, -1, x) - 6H{-r, -r, 0, -1, x) 
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+3i/(-r, -4, -r, -1, x) + 2H{-r, -4, -4, -r, x) 
-3H{0, -r, -r, -1, - 2H{0, -r, -4, -r, x) 

+2H{0,0,-r,-r,x)Y (209) 
The above MI has a simple UV pole coming from the sub-divergence in the bubble. 




^ J V^VrVs 



(210) 
(211) 



where: 



aF, 



(17) 



1 r 

X 



C(2)i7(0, l,a;)-2i7(0, 1,0,-1, x) + 2//(0,r,r, 0,x) . (212) 



4.4 Topology t = Q 




(3/ / 



(213) 
(214) 



where: 



2 77(18) 







C(2) (^i/(0, 1, x) + H{1, 1, x) + H{1, c, x) + H{l,c, x) 

+H{0, c, x) + H{0, c, x)^ - 2H{1, 1, 0, -1, x) - 2/7(0, 1, 0, -1, x) 
-3H{0,c,0,-l,x) -3H{0,c,0,-l,x) -2H{l,r, r,0,x) 
-2H{0,r,r,Q,x) -3H{l,c,Q,-l,x) -3H{l,c,Q,-l,x) 



X 



-i-H{r, 0, 1) H{0, c, x) - H{0, c, x) + H{1, c, x) 



-H{l,c,x) 



(215) 



The index r appears in the GHPLs only or 2 times, so the coefficients of the related 
terms do not contain radicals. 




Oj / 



(216) 
(217) 
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where: 



a'Fo^''^ = ^{m[H{0,-l,x)-H{0,l,x)-H{l,l,x) + H{l,-l,x) 
-2H{0, -1, 0, -1, x) + H{0, -1, 0, 0, x) + 2H{0, 1,0, -l,x) 
-2H{1, -1, 0, -1, x) + H{1, -1, 0, 0, x) + 2H{1, 1, 0, -1, x)| . (218) 

Because of analogous considerations to the previous ones, the above MI is expressed 
in terms of ordinary HPLs. 




where: 



2(4-D) f 'J^Dj^^Dj. 1 



(219) 
(220) 



2 17.(20) 







: 12C(2)iJ(-r, -1, x) - 6H{-r, -1, 0, -1, x) 



x^ x{x + 4) t 
+6i/(-r, -1, 0, 0, x) - 12H{-r, -r, -r, -1, x) 
+8H{—r, —r, —r, 0, x) + 8H{—r, —r, 0, — r, x) 

+4:H{-r, 0, -r, -r, x) - 2H{-r, 0, 0, -1, x) 



(221) 



where: 



(Jj / 



1 



(222) 
(223) 



2 7r(21) 







x^ 



6C(2)iy(l, 1, x) + 6H{1, r, r, 0, x) - 2H{1, 0, 0, -1, x) 



+if(l,0,l,0,x) - 12i7(l, 1,0,-1, x) + 4iy(l, 1,0,0, x) . (224) 

The masses of the vector bosons exchanged in the t channel completely cut-off the 
infrared singularities, so the above MI is IR (as well as UV) finite. As noted in our 
previous work P], a non-zero mass on the outer boson line is already sufficient to 
completely screen the IR singularities. 

The above MI has been computed in [21j by fitting a small momentum expan- 
sion to an assumed form for the exact expression. We have compared the first 15 
terms of the small momentum expansion in with an analogous expansion of 
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our expression, finding complete agreement. The first few terms of tlie small mo- 
mentum expansion of are given in the next section. In [221 ^ leading-twist 
large-momentum expansion of the above MI has been presented, based on the sep- 
aration of the loop space in leading IR regions and a related approximation on the 
integrand. The above result is in agreement with a preliminary large momentum 
expansion of our expression [THj . 

5 Reducible six-denominator amplitudes 




a J 



(225) 
(226) 



where: 



C(2)/f(0, 1, x) + C(2)i/(1, 1, x) - HiO, 1, 0, 0, x) 



+i/(0, 1, 1, 0, x) - if (1, 1, 0, 0, x) + i/(l, 1, 1, 0, x) 



(227) 



In the above amplitude, the photon is emitted by the W boson internally to the 
triangle. It is therefore "trapped" and cannot propagate for large distances. The 
consequence is that the above amplitude is IR finite, as expected on the basis of 
physical intuition. 




1 



^ ^ i=—l 



where: 



2 p(23) 



2H{0, 0, -1, x) + H{0, 1, 0, x) + 2H(1, 0, -1, x) 







+if(l,l,0,x 

-- k(2)(if(0, 1, x) + H{1, 1, x)) - 8H{0, 0, -1, -1, x) 
X L 

+H{0, 0, 1, 0, x) - H{0, 1, 0, 0, x) + 4:H{0, 1, 1, 0, x) 
-8/7(1, 0, -1, -1, x) + H{1, 0, 1, 0, x) - H{1, 1, 0, 0, x) 

+4if(l,l,l,0,x) . 



(228) 
(229) 



(230) 



(231) 
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The above amphtude corresponds to a photon emitted hy a W boson outside the 
triangle. In this case, the photon can propagate to large distances and a coUinear 
singularity is generated, characterized by the presence of the simple 1/e pole. As 
anticipated in the introduction, a qualitative analysis of the infrared singularities 
can be done by shrinking all the internal massive lines to a point: in this limit the 
diagram factorizes into a massless bubble evaluated at the light-cone momentum 
P2 times a massless bubble evaluated at the general momentum q. The coUinear 
singularity originate from the former bubble, representing the evolution of a jet 
formed by an initial particle. 




^ ^ i=-2 



(232) 
(233) 



where: 






1 




X 




1 - 


X - 




1 - 


X - 



l + H{0,x)-H{-l,x) 



1 



(234) 

H(-l,x) + H(0,-l,x) , (235) 



C(2) - 2H{0, x) + 2H{-1, x) - H{0, 0, x) + 4H{-1, -1, x) 
1 r 



-4:H{0,-l,x) 



2H{-1, x) + 4:H{-1, -1, x) + 2/7(0, -1, x) 



+2H{r, r, 0, x) + 4H{0, -1, -1, x) + 2H{0, 0, -1, x) 

A- X 



x^/x{4: — x) 



H{r,0,x) 



(236) 



4 + 2C(2) - 2C(3) - 4//(0, x) + C(2)//(0, x) + 2C(2)//(1, x) 

+4:H{-1, x) + 2C(2)/7(-l, x) - 2H{0, 0, x) - 8H{0, -l,x) 

+8H{-1, -1, x) - H{0, 0, 0, x) + 16H{-1, -1, -1, x) 

-6H(-1, 0, -1, x) - 16H(0, -1, -1, x) + 2H(0, 0, -1, x) 

1 



-AH{1, 0, -l,x) + 3H{r, r, 0, x) 



+ 



x^ 



AHi-l,x) + 2Ci2)Hi-l,x) 



-2C(2)i7(l, a;) - 2C(2)ff(0, 1, x) + 4iJ(0, -1, x) + 2C(2)/7(0, -1, ,t) 
-2C{2)H{r, r,,T) + 8i/(-l,-Lx) + 4i7(0, 0,-1,, t) + 4iJ(l, 0,-1, x) 
+ 16if(-l, -1, -1, x) - 6i/(-l, 0, -1, x) + 8i7(0, -1, -1, x) 
+2^r(r, r, 0, 0, x) + 2iy(r, 0, r, 0, x) - 4H{r, 4, r, 0, x) 
+16iJ(0, -1, -1, -l,x) - 6H{0, -1, 0, -1, x) + 8H{0, 0, -1, -1, x) 

+4:H{0, 0, 0, -1, x) + 4:H{0, 1, 0, -1, x)^ ~ ^ 



a;iya;(4 — x) 
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-2H{r, 0, x) - H{r, 0, 0, x) - H{0, r, 0, x) + 2H{4, r, 0, x) 



(237) 



The above amplitude has a double pole in e coming from the following region: in 
the inner triangle a large momentum k2 flows while in the external triangle a soft 
momentum ki flows, or /cf -C In this limit, the diagram factorizes into a massless 
triangle with the well-known double IR pole (soft x collinear), times a coefficient 
functions given by the inner triangle, which is effectively point-like. 



6 Small momentum expansion of six denominator 
amplitudes 

In this Section we present the small momentum expansions \s\ -C of all the 
6-denominator diagrams, i.e. expansions in powers of x and L — logx up to first 
order in x included. The expansion of the GHPLs for a small value of the argument 
-C 1 is obtained in the following way: 

• we explicitly write the GHPL as a repeated integration over the basis functions, 
as for example: 



H{l,r,r,0;x)= / / 
Jo l-xi Jo 



dxo 



rx2 

Jo 



dx^ 



•>/x2(4 - X2) Jo a/x3(4 - X3) 



logxs; 



(238) 



we expand the basis functions in powers of x up to the required order in x^^, 
such as for instance: 



7^(4 - x) 2Vi 16 



(239) 



we integrate term by term the expanded functions. This involves in general 
the integration of functions of the form x"^ loga;*^, with q integer or half integer 
and k integer. 



The expansions are given below. 



where: 




2\ 2e 2 

5^a;M; + C»(x3), 

j=0 



a'A'o 



= 4 + V3//(r,0,l) -L, 



i^ + lc(2) + ^i/(r,0,l)-l^L, 
9 2^ ^ 4 ^ ' ' ^ 24 

12509 2,, , \/3,,, ^ , 197 

+ -C 2 + —H(r, 0, 1) . 

16200 3^^ ^ 2 ^ ' ' ^ 540 



17. 



The factors 1/x and l/y/x are clearly not expanded. 



(240) 

(241) 

(242) 

(243) 
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where: 




j=0 



(244) 



^2 rO 



3-2L + -L^, 

8 3^^ ^ 108 36 



(245) 
(246) 
(247) 



where: 




2\ 2c 2 
j=0 



a;2 



(248) 



1 + 2C(2)-L+-L2, 

-^-C(2) + ^L-iL^ 
24 ^ 12 4 ' 



2 — 



827 1 
2160 ^ 2^^^'' 



61 
360 



L + -L\ 



(249) 
(250) 
(251) 



where: 




(252) 



/,2p0 

,2 TTiO 



-4 + 3C(2) -L + L\ 



.HH! + Hc(2) 

600 4 ^ 120 



L H 

12 



(253) 

(254) 

(255) 



where: 




(256) 



3 - C(2) - 2L + -L^ 

^ _ ^/•('2) - -L + -L^ 
16 4^^ ^ 4 8 ' 

251 11 , , 107, 11 2 

C 2 L + —L^ 

216 18^^ ' 108 36 



(257) 
(258) 
(259) 
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where: 



where: 




/ 7=0 i=-l 



-6 + AC{2) + 3L-h\ 

-4 + 4^' 

11 . . 21 3 p 

12 18 

216 9 ^ 9 36 




(260) 

(261) 

(262) 

(263) 
(264) 
(265) 
(266) 



^ ^ 7=-l i=-2 





1, 




(268) 




1-L, 




(269) 




l-c(2)-L+lL^ 




(270) 




1-C(2)-2C(3)+C(2)L-L + 


2 6 


(271) 


aVo"^ = 


1 

4' 




(272) 




35 1 , 

^ L, 

72 12 ' 




(273) 


aV(J = 


_35^13 ^L- 

432 12^^ ^ 24 24 ' 




(274) 


aVf' = 


1 

^18' 




(275) 


aV° = 


209 1 , 

\ L, 

675 180 ' 




(276) 


aV^ = 


141091 1 , , 23 ^ 
162000 180^^ ^ 1080 


360 ' 


(277) 




1 

48 ' 




(278) 


aV° = 


11063 1 ^ 
78400 ' 1680 ' 




(279) 




105455939 47 17 , 
296352000 560^' ' 6720 


- 1 
3360 


(280) 
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7 Large momentum expansion of six denominator 
amplitudes 



In this Section we give the asymptotic expansions for \s\ S> rri^ of all the 6- 
denominator scalar integrals, i.e. the expansion in powers oi 1/x and L — \ogx 
up to the order included. These results are relevant for the study of the struc- 
ture of the infrared logarithms coming from multiple emission. 



(281) 



where: 



-8a, + ^C'(2) + 2C(2) V 2 - ^ log' 2 , 

-3 - 2C(2) + C(3) - 2C(2)L -3L- - , 
91 1 1 1 1 



(282) 
(283) 
(284) 
(285) 



and where a, — Li4(l/2). 



where: 




2\ 2e 4 



E- 



O 



(286) 



a'G\ = 
a^G\ = 



-C'(2) , 
10^ ^ ^ ' 

-2-C(2)-C(3)-LC(2)-2L-L2 



^^C(2)4c(3)4c(2)L+^L + 1l^-1l3^ 



(287) 
(288) 
(289) 



+ 1^(2) - 1c(3) - ic(2)L + —L + - . (290) 
162 18^^ ' 3^^ ' 2/^ ^ 108 18 9 ^ ^ 



where: 




2\ 2e 4 2 

j=2 i=-2 



+ 



X'^ 



(291) 



-2C(3) , 



(292) 
(293) 
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(294) 




-4 + 2C(2)-4C(3) + 2C(2)L-4L 



2L2 




(295) 






y-^C(2)-2C(3) + C(2)L+^L 




(296) 



(297) 




13 11. 1 



(298) 



i 

108 18 6 





(299) 
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Conclusions 



We have presented the exact analytic evaluation of the 25 master integrals con- 
taining 2 and 3 massive propagators entering the planar amplitudes of the 2-loop 
electro weak form factor. While the reduction to master integrals does not present 
any new element with respect to our previous computation ^ and is done with the 
same algorithm, the analytic evaluation of the master integrals requires a non-trivial 
extension of the harmonic polylogarithm theory. The presence of 2 massive parti- 
cles in the s or in the t channel opens indeed thresholds and pseudothresholds in 
s = ±4m^ respectively, in addition to the old ones in s = 0, ±m^. 

The generalization of the 1-dimensional harmonic polylogarithms has basically 
required: 

• the introduction of new basis functions, in addition to the usual one, involving 
complex constants and radicals; 

• a set of recursion relations to take the integrals with semi-integer powers com- 
ing from the evaluation of the master integrals to a unique form fixed by the 
basis function choice. 

The basic properties of the ordinary harmonic polylogarithms are maintained by the 
generalization. 

The small momentum expansion of all the 6-denominator amplitudes has been 
obtained by means of a series expansion of the basis functions. 

We could also obtain the large momentum expansion of all the six-denominator 
amplitude involving only ordinary harmonic polylogarithms. 

We compared our results with those present in the literature usually in the form 
of resummed small momentum expansions or truncated large momentum expansions, 
finding complete agreement. 

In order to complete the evaluation of the master integrals, 3 steps are still to 
be taken: 

• as explained in Sectional the transformation x — > 1/x requires the knowledge 
of all the H{w; a;)'s in x = 1. The H{w; l)'s have to be expressed in terms of 
a minimal set of transcendental constants |18j : 
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• the evaluation of the master integrals related to the crossed ladder topology. 
The reduction, using both the numerical-indices method and the symbolic 
method, shows that this topology has 3 master integrals. The resulting system 
of 3 differential equations cannot be completely triangularized by means of the 
techniques discussed in this work, but can be split into a second-order and a 
first-order differential equations [TH] . 

• The numerical evaluation of the generalized harmonic polylogarithms, which 
does not seem to have specific difficulties with respect to the ordinary case. 
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A One-loop master integrals 

In this Appendix we present the results for the 1-loop master integrals containing 2 
massive propagators. We have recomputed them with the method of the differential 
equations described in the main body of the paper in terms of GHPLs. In the 
case of the bubble, we found that it was necessary to push the e expansion up to 
third order included. In our previous work P we gave the expressions of the 1-loop 
master integrals containing at most 1 massive propagator. The above amplitudes 
are necessary for the computation of the factorized 2-loop master integrals, which 
are presented in the next section. 



A.l Bubble 





(300) 



where: 



Bo 



2 - 



1 



y/x{x + 4) 



X + 4 



H{-r,x) , 



(301) 
(302) 
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= 4 - 


B2 


= 8 

A 


B3 


= 16- 



X + 4 



^/x{x + 4) 
X + 4 



2H{-r,x) -H{-4:,-r,x) 



(303) 



X + 4 



a/x(x + 4) 
-if (-4, -4, -4, -r,x 



AH{-r,x)-2H{-A~r,x) + H{-A-A-r,x) , (304) 
if (-r, x) - 4ff (-4, -r, x) + 2ff (-4, -4, -r, x) 

(305) 



A. 2 Vertex 




(4-D) / ^D^ 



[(pi - kf + a] [(p2 + kf + a] 
(^)'i:^-V, + 0(e'). 

^ ^ i=0 



(306) 



where: 



af^o 
aKi 
aK2 



—H(—r, —r, x) , 

X 

9 r 

H{—r, —4, — r, x) — if (0, — r, — r, x) 
H{—r, —A, —A, —r, x) — if (0, — r, —4, — r, x) 



2 

X L 



(307) 
(308) 

(309) 



+if(0,0, -r, -r,x) 

B Factorized master integrals 

In this Appendix we give the expressions of the factorized 2-loop master integrals, 
i.e. of the Mis in which the 2 loops do not have common propagators. One has only 
to multiply the 1-loop master integrals representing the separated subdiagrams and 
convert products of GHPLs ff (a; x)H{b; x) into linear combinations of ff 's by using 
the algebra-identity in [2]. 



B.l Topology t = 2 



+ (311) 
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where: 



.(25) 



-2 



^(25) 



.(25) 



(25) 



.(25) 



1, 

2, 
3, 
4, 
5. 



(312) 
(313) 
(314) 
(315) 
(316) 



The above amphtude appears, in general, in the reduction of all the amplitudes 
having at least 1 massive propagator in anyone of the 2 loops. 

B.2 Topology t = 4 



= (^)"e^'';'"'+o(^=). 



where: 

^(2^6) 
^(26) 



1, 

A-H(0;x) 



V^(4 + 



X] 



--H{-r;x), 



(317) 
(318) 



(319) 
(320) 



.(26) 



= 12-C(2)-4//(0;x)+//(0,0;x)- 



+i7(-4, -r; x) + H{Q, -r; x) 



(4 + a;) 
^/x{ITx) 



-4:H{-r;x)+H{-r, 0;x) 
(321) 



(26) 



-12H{-r; x) +C{2)H{-r; x) +AH{-r, 0; x) - H{-r, 0, 0; x) 



.(26) 



32 - 4C(2) - 2C(3) - 12H{0; x) + H{0; x)C{2) +AH{0, 0; x) - H{0, 0, 0; x) 
I (4 + ^) 
V^a;(4 + a;) 

+4ii'(-4, -r; x) - ^^(-4, -r, 0; x) - ^^(-4, -4, -r; x) - ^^(-4, 0, -r; x) 
+4//(0, -r; x) - H{0, -r, 0; x) - H{0, -4, -r; x) - H{0, 0, -r; x)] , (322) 

80 - 12C(2) - 1(2(2) _8^(3) - 32iJ(0; x) +4C(2)i7(0; x) +2C(3)i/(0; x) 
+12ii"(0, 0; x) - C(2)i/(0, 0; x) - 4:H{Q, 0, 0; x) + i/(0, 0, 0, 0; x) 
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(4 + a;) 



-32H{-r; x)+2(2C(2) +C(3))//(-r; x) + 12H{-r, 0; x) 



a/x(4 + x) 

-C{2)H{-r, 0; - 4/J(-r, 0, 0; + H{-r, 0, 0, 0; x) + 12i/(-4, -r; x) 
-C{2)H{-A, -r; - 4iJ(-4, -r, 0; x) + if(-4, -r, 0, 0; x) 
-4i7(-4, -4, -r; x) + if(-4, -4, -r, 0; x) + H{-A, -4, -4, -r; x) 
+H{-A, -A, 0, -r; x) - 4H{-A, 0, -r; x) + H{-A, 0, -r, 0; x) 
+H{-A, 0, -4, -r; x) + i/(-4, 0, 0, -r; x) + 12if(0, -r; x) 
-C{2)H{0, -r; x) - 4i7(0, -r, 0; x) + H{0, -r, 0, 0; x) - 4i/(0, -4, -r; x) 
+H{0, -4, -r, 0; x) + H{0, -4, -4, -r; x) + H{0, -4, 0, -r; x) 
-4//(0, 0, -r; x) + H{0, 0, -r, 0; x) + H{0, 0, -4, -r; x) 



+1^(0, 0, 0, -r; x) 



(323) 



where: 



^(27) 



(27) 




^ ' i=— 1 



1^121^13 



H{l,0,x), 

C{2)H{1, x) + H{1, 0, x) + i/(0, 1, 0, x) - i/(l, 0, 0, x) 



+//(!, l,0,x) 
1 r 



(324) 
(325) 



(326) 



(327) 



X 



(C(2) + 2C(3))iJ(l, x) + (l - C(2))i/(1, 0, x) + C(2)if(l, 1, x) 

+C(2)if(0, 1, x) + H{1, 1, 0, x) + H{0, 1, 0, x) - H{1., 0, 0, x) 
+H{0, 0, 1, 0, x) - H{0, 1, 0, 0, x) + H{0, 1, 1, 0, x) + H{1, 0, 0, 0, x) 



+H{1, 0, 1, 0, x) - H{1, 1, 0, 0, x) + H{1, 1, 1, 0, x) 



(328) 



B.3 Topology t = 5 




^ ' i=— 1 



10^^15^^16 



(329) 
(330) 
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where: 



(28) 



aF, 



(28) 



H{—r, —r, x) 



(331) 



2H(—r, —r, x) — H{—r, —r, 0, x) — H(—r, —4, — r, x) 



-H{—r, 0, — r, x) 



(332) 



(28) 



(4 - C(2))if(-r, -r, x) - 2H{-r, -r, 0, x) - 2H{-r, -4, -r, x) 

-2H{-r, 0, -r, x) + ^r(-r, -r, 0, 0, x) + H{-r, -4, -r, 0, x) 
+H{-r, -4, -4, -r, x) + H{-r, -4, 0, -r, x) + H{-r, 0, -r, 0, x) 



-H{-r, 0, -4, -r, x) + //(-r, 0, 0, -r, x) 



C Reducible two-loop amplitudes 



(333) 



In this Appendix we present the expressions of some interesting 2-loop amphtudes 
which can be reduced to the Mis given in the present paper and in pp. 

C.l Topology t = 3 

This amplitude reduces to the product of tadpoles coming from the contraction of 
its massless line. 



where: 



i=-2 



.(29) 



a 

p(29) 
^0 



a 



(29) 



.(29) 



-1, 

-3, 

-r, 

-15, 
-31. 



(334) 
(335) 



(336) 
(337) 
(338) 
(339) 
(340) 
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C.2 Topology t = 5 



(30) 




2(4-D) 



1 



i=-l 



where: 






-1/ 




a; 




1 - 


X . 



(341) 

(342) 

(343) 



2C{2)H{1, x) + H{0, -1, x) + H{0, 1, 0, x) + H{1, 0, x) 

-2H{l.,0,-l,x) - H{l,0,0,x) + H{l,l,0,x) 
1 



+ 



C(2)-2//(0,-l,x) 



(l-x) 

^ ^C(2)i/(1, x) + C(3)if(l, x) + 3H{0, -1, x) + H{1, 0, x) 



(344) 



X 



-C(2)i7(l, 0, ,t) + a2)H{l, 1, x) + 2C(2)i7(0, l,x) + F(0. 0, -1, x) 
+H{0, 1, 0, x) - H{1, 0, 0, ,t) + H{1, 1, 0, x) - AH{0, -1, -1, x) 
+i7(0, 0, 1, 0, x) - 2H{0, 1, 0, -1, x) - H{0, 1,0,0, x) 
+H{Q, 1, 1, 0, x) + H{1, 0, 0, 0, x) + H{1, 0, 1, 0, x) 

-H{1, 1, 0, 0, x) + H{1, 1, 1, 0, x) + 8H{1, 0, -1, -1, x) 
1 



+ 



(l-x) 



3C(2) - C(3) - 6H{0, -1, x) + 8i/(0, -1, -1, x) . (345) 




2(4-D) 



1 



^ ^ i=— 1 



(346) 
(347) 



where: 












--H{l,0,x), 

X 




(348) 




1 


C{2)H{l,x) + 2H{l,0,x)- 


H{l,0,0,x) + H{l,l,0,x) 


, (349) 




1 

X 


'2a2)H{l,x) + 2C{3)H{l,x 


) + 4H{l,0,x) + am{hl,x) 



-2C(2)//(0, 1, x) - C(2)//(l, 0, x) + 2H{1, 1, 0, x) - 2H{1, 0, 0, x) 
-2H{0, 1, 1, 0, x) + H{1, 0, 0, 0, x) + H{1, 0, 1, 0, x) 

-H{1, 1, 0, 0, x) + H{1, 1, 1, 0, x)] . (350) 
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The above 2 amplitudes have a simple ultraviolet pole coming from the nested 
bubble, containing 1 massive denominator. The only difference is that the bubble is 
inserted in the s-channel in the former diagram and in the t-channel in the latter. 
The coefficient of the simple pole is indeed the same. 




1 



^ ' i=—3 



(351) 

(352) 



where: 



aF 



(32) 



1 

X 



aF 



(32) 
-2 



-H{Q,x), 

X 

Hc(2)-//(0,0,x) 



aF, 



(32) 



X 

-V 

X . 



(353) 
(354) 
(355) 



S- 2C(3) +H{r, r, 0, x) - 4H{0, x) + C{2)H{0, x) - H{0, 0, 0, x 

2(4 - x) 



aFl 



(32) 



X^yx{4: — x) 
1 

X 



H{r,0,x), 



9 



(356) 



16 + 4C(2) + ^C'(2) - 2C(3)//(0, x) - (4 - C(2))//(0, 0, x) 
-C{2)H{r, r, x) + 4H{r, r, 0, x) - H{0, 0, 0, 0, x) + H{r, r, 0, 0, x) 
+H{r, 0, r, 0, x) - 2H{r, 4, r, 0, x) - H{0, r, r, 0, x) 
2(4 - x) 



4 

a;(4 — x) 

+2i7(4,r,0,a;) 



C{2)H{r, x) - H{r, 0, 0, x) - H{0, r, 0, x) 



(357) 



In the IR limit, the above amplitude factorizes into the product of a massless 1-loop 
triangle times a vacuum bubble with 2 masses. The triple pole is the product of 
a double IR pole coming from the triangle and a simple UV pole coming from the 
bubble. 




^2(4-D) J ^Dj^^^Dj^^^ 



1 



^ ^ i=—l 



(358) 
(359) 
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where: 



aFr 



(33) 



— if (0,0,-1, x), 

X 

^ \({2)H{0, 1, x) - H{0, 0, 0, -1, x) - 4:H{0, 0, -1, -1, x) 



X 



-2H{0,l,0,-l,x) 



(360) 



(361) 



According to IR power counting, one can shrink all the massive lines in the ampli- 
tude. This shows that the above amplitude has a simple collinear pole associated 
to the evolution of the lower external line. 



D Scalar diagrams of self-energy type 

In this Appendix we present the results for the scalar diagrams of self-energy in- 
sertion type (see Fig. |2I). As explained in detail in Section |21 these diagrams are 
effectively 5- and 4-denominator amplitudes and are all reducible to the Mis by 
means of the ibps identities. Their analytic expressions are given for completeness. 




1 



i=-2 



(362) 
(363) 



where: 



2 r(34) 



a'F 



a i^o 



1 -if(0,x) 



(1-x) 



if(0,x) 



(364) 
(365) 



2x^ X 



- - C(2) - if(0, x) - H{-l,x) + H{0, 0, x) 



+ ^H{0,-l,x)-^H{l,0,x) 



-H{0,0,x) + H{l,0,x) 



'l-x) 



a2) + H{0,x)+H{-l,x) 



(l-x) 



1 - 



(l-x). 



:C(2) 



-if(0,-l,x) 



(366) 



x^ 



7 



H{-l,x) -2H{-l,-l,x) 



11 

T 



C(2)-2C(3)-if(0,x) 



+C(2)if(0, x) - ^-H{-l,x) - ^amih x) + if (0, 0, x) - ^ii(l, 0, x) 
+4ii(-l, -1, x) - ^ii(0, -1, x) - ii(0, 0, 0, x) - 2ii(0, -1, -1, x) 
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+^H{0, 0, -1, x) - ^H{0, 1, 0, x) + 2H{1, 0, -1, x) + ^H{1, 0, 0, x) 



--H{l,l,Q,x) 



1 



(l-x) 



C(2)-2C(3)-(l-C(2))i/(0,a;) 



-C(2)i/(1, x) - x) + i/(0, 0, x) - H{1, 0, x) + 4i7(-l, -1, 

-4i7(0, -1, x) - i/(0, 0, 0, x) - H{0, 1, 0, x) + H{1, 0, 0, x) 



-H{l,l,0,x) 



+ 



(l-x) 



1 



(l-x) J 



^C(2) - ^C(3) 



--i/(0,-l,x) + 4i^(0,-l,-l,x) 



(367) 



x^ 



35 

-i/(-l,x) + C(2)i^(-l,x)-7i/(-l-l,x) + 8i/(-l-l-l,x) 



-3ii'(-l,0,-l,x) 



1 

+ - 

X 



43 Q 

-- + -C(2)2 + 2C(3) + (1-C(2) 



-2C(3))//(0, x) + -jH{-l, x) + 2C(2)i/(-l, x) + -C(2)//(l, x) 

+2C(3)//(1, x) - (1 - C(2))i/(0, 0, x) - UH{-1, -1, x) 

49 5 15 

+ -i7(0, -1, x) - C(2)i/(0, -1, x) + -C(2)i/(0, 1, x) + -i/(l, 0, x) 

-^C(2)i/(1, 0, x) + ^C(2)i/(1, 1, x) + H{0, 0, 0, x) - 5H{0, -1, -1, x) 

+^ii'(0, 0, -1, x) + ^i/(0, 1, 0, x) - H{1, 0, -1, x) - ^i/(l, 0, 0, x) 

7 

+16//(-l, -1, -1, x) - 6i/(-l, 0, -1, x) + -//(1, 1, 0, x) 
-i7(0, 0, 0, 0, x) + 2H{0, 0, -1, -1, x) - 8H{0, -1, -1, -1, x) 
+3H{Q, -1, 0, -1, x) - ^H{0, 0, 0, -1, x) + ^H{Q, 0, 1, 0, x) 

-2H{0, 1, 0, -1, x) - ^i/(0, 1, 0, 0, x) + ^H{0, 1, 1, 0, x) 

+8H{1, 0, -1, -1, x) + ^H{1, 0, 0, 0, x) + ^H{1, 0, 1, 0, x) 



-^i/(l,l,0,0,x) + ^i7(l,l,l,0,x 



1 



6C(2) - ^C^(2) 



(l-x)L 

-4C(3) - HiO, x) + C(2)i/(0, x) + 2C(3)i^(0, x) - 2C(3)i/(l, x) 
35 

--H{-1, x) - 2C(2)//(-l, x) + H{0, 0, x) - C(2)//(0, 0, x) 

-//(I, 0, x) + C(2)i/(1, 0, x) - C(2)i/(0, 1, x) - C(2)//(l, 1, x) 
+14//(-l, -1, x) - UH{0, -1, x) - H{0, 0, 0, x) - H{0, 1, 0, x) 
-2H{1, 0, -1, x) + H{1, 0, 0, x) - H{1, 1, 0, x) + 16//(0, -1, -1, x) 
-2H{0, 0, -1, x) - 16H{-1, -1, -1, x) + 6H{-1, 0, -1, x) 
+i/(0, 0, 0, 0, x) - H{0, 0, 1, 0, x) + H{0, 1, 0, 0, x) 
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-H{0, 1, 1, 0, x) - H{1, 0, 0, 0, x) - H{1, 0, 1, 0, x) 
+H(l,l,0,0,x)- H(l,l,l,0,x) 



(1-x) 
35 



'^5 Q 7 

-C(2) + -C^(2)--C(3)+C(2)i/(0,l,a;) 



-— if(0,-Lx) - 2C(2)H{0,-l,x) + 14i/(0,-l,-l,,T) 
-16/^(0, -1, -1, -1, x) + 6H{Q, -1, 0, -1, x) - 2H{Q, 0, 0, -1, x) 

-2H{0,l,0,-l,x)\ . (368) 



where: 



2 c.(35) 



-2 



2 r(35) 



2 E.(35) 







2 c.(35) 




^ J VlVJ)^ 



(369) 
(370) 



2 1 

6x ' 

^ [l + if(0,x 



1 



x^ 



g^i/(0,x), 
^ - 2C(2) + ^H{0, x) + 2H{0, 0, x) 



(371) 
(372) 



1 

+ - 

X 



~H{0,x)-^H{0,0,x) 



3^/x{4 



X 



x^ X ~'~ 2 



H{r, 0,x)(373) 



1 



x^ 



^ + ^C(2) - 4C(3) - ^H{0, x) + 2C(2)i^(0, x) - ^H{0, 0, x) 



-2ii"(0,0,0,x) 
4 



^.>)-lc(3)-(i->,.(0.) 



gif(0, 0, x) - ^H{0, 0, 0, x) + ^H{r, r, 0, x 



1 



_4 1 

~x ^ 2 



3yx(4 - x) 

C{2)H{r, x) - H{r, 0, 0, x) - H{0, r, 0, x) + 2if(4, r, 0, x 



x^ 



+ - 



x^ 



4 



i/(r,0,x)L 



(374) 



"20 _ 1 _ ^ 

x^ X 

58 \ 2 

- - 2C(2) j (0, 0, x) - -H{0, 0, 0, x) + AH{r, r, 0, x) 
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1 

+ - 

X 



968 26^,^, 3 ,,^,2 8. /ox 



-2//(0,0,0,0,x) 

-^amir, r, x) - ^H{0, 0, 0, x) + ^H{r, r, 0, x) - ^H{0, 0, 0, 0, x) 
+^H{r, r, 0, 0, x) + ^H{r, 0, r, 0, x) - ^H{0, r, r, 0, x) 



-H{r,4:,r,0,x) 



X' 



X] 



^Cm{r:X) + ^mH{r,x) 



128 8 8 16 

+— i/(r, 0, x) - -a2)H{r, 0, x) - -C(2)i/(0, r, x) + y C(2)i/(4, r, x) 

-^//(r, 0, 0, x) - ^H{0, r, 0, x) + ^^^(4, r, 0, x) + SH{r, r, r, 0, x) 

+ ^H{r, 0, 0, 0, x) + ^H{0, r, 0, 0, x) + ^H{0, 0, r, 0, x) 

16 16 16 

-y i/(0, 4, r, 0, x) - -HiA, r, 0, 0, x) - -H{A, 0, r, 0, x) 



32 

+-i/(4,4,r,0,x) 



1 

+ - 



1 



X ^yx{A - x) 



^a2)H{r,x) + ^Cmir,x) 



+ ^H{r, 0, x) - ^mH{r, 0, x) - ^C(2)i^(0, r, x) + ^C(2)i^(4, r, x) 

-^H{r, 0, 0, x) - ^H{0, r, 0, x) + ^iJ(4, r, 0, + 4/f(r, r, r, 0, x) 
9 9 9 

+^H{r, 0, 0, 0, x) + ^H{0, r, 0, 0, x) + ^H{0, 0, r, 0, x) 
-^HiO, 4, r, 0, x) - ^H{4, r, 0, 0, x) - ^H{4, 0, r, 0, x) 



+yi/(4,4,r,0,x) 



+ 



X] 



^am{r,x)-\cm{r,x) 



-^H{r, 0, x) + ^amir, 0, x) + ^C(2)//(0, r, x) - ^C(2)//(4, r, x) 
-^i/(r, 0, 0, x) - ^H{0, r, 0, x) + ^if(4, r, 0, - ^i/(r, r, r, 0, 
-^//(r, 0, 0, 0, x) - ^H{0, r, 0, 0, x) - ^H{0, 0, r, 0, x) 
+ Ih{Q, 4, r, 0, x) + ^H{A, r, 0, 0, x) + ^H{A, 0, r, 0, a;) 



--//(4,4,r,0,x) 



(375) 



1 



(376) 
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where: 



(l-x) 



H{0,x) 



X (1 — x) 



+H{l,0,x 
1 



C{2) + 2H{0,x) -2H{0,0,x) 



(377) 



(378) 



(379) 



X 



3C{2)H{1, x) + 6H{1, 0, x) + 3H{0, 1, 0, x) - QH{1, 0, 0, x) 
1 

AH{0, 0, x) + 2/7(1, 0, x) + AH{0, 0, 0, x) + H{0, 1, 0, x) 



+3/7(1, l,0,x) 



2C(2) + 3C(3) + 4/7(0, x) + C(2)/7(l, x) 



-2/7(1, 0,0,x) + 77(1, 1,0, x) , (380) 

-- [6C(2)//(1, x) + 9C(3)//(1, x) + 12//(1, 0, x) + 3C(2)//(0, 1, x) 
X L 

+3C(2)//(1, 1, x) + 6//(0, 1, 0, x) + 6//(l, 1, 0, x) - 12//(1, 0, 0, x) 
+3//(0, 0, 1, 0, x) - 6//(0, 1, 0, 0, x) + 3//(0, 1, 1, 0, x) 
+12/7(1, 0, 0, 0, x) + 37/(1, 0, 1, 0, x) - 6//(l, 1, 0, 0, x) 

1 

-8C(3)//(0, x) + 2C(2)//(1, x) + 3C(3)//(1, x) - 8//(0, 0, x) 
+4//(l, 0, x) + C(2)//(l, 1, x) + C(2)//(0, 1, x) + 8//(0, 0, 0, x) 
+2//(0, 1, 0, x) - 4//(l, 0, 0, x) + 2//(l, 1, 0, x) - 8//(0, 0, 0, 0, x 
+//(0, 0, 1, 0, x) - 2//(0, 1, 0, 0, x) + //(0, 1, 1, 0, x) 
+4//(l, 0, 0, 0, x) + 7^(1, 0, 1, 0, x) - 2//(l, 1, 0, 0, x) 

+//(l,l,l,0,x) 



+3//(l,l,l,0,x) 



4C(2) + ^C'(2) + 6C(3) + 8//(0,x) 



(381) 



The following diagram is reducible to a combination of two different 5-denominator 
diagrams, as explained in Eqs. (5,6) of JLJ, by simple partial fractioning. 



(382) 
(383) 




2(4-D) f j^i^^ 2)^fc2— — — ^ 



12^^13/^14 



where: 



X 



(384) 



62 



—H{0,x), 

X 

nC(2)-i/(0,0,a;)-i/(l,0,a;) 



X 



(385) 

(386) 



8 - 2C(3) - (4 - C(2)) //(O, x) + (C(2) + V?,H{r, 0; 1)) (1, x) 

+2//(l, 0, x) - H{0, 0, 0, x) + H{0, 1, 0, x) - H{1, 0, 0, x) 

1 2iA — x\ 

+H{1, 1, 0, x) - 3H{1 + r, r, 0, x) + -^==L=H{r, 0, x) , 

J x(4 — X) 

-i|l6 + 4C(2) + ^C'(2) - 2C(3)i/(0,a;) - 2(c(2) + C(3))i^(l, x) 

- (4-C(2)) (//(O, 0, x)+H{l, 0, x)) -C(2)//(0, 1, x) -C(2)//(l, 1, x) 

+3C(2)i7(l + r, r, x) - 2H(0, 1, 0, x) + 2H(1, 0, 0, x) - 2H(1, 1, 0, x) 
+6//(l + r, r, 0, x) + 6H{r, r, 0, x) - 3H{1 + r, r, 0, 0, x) 
-3H{1 + r, 0, r, 0, x) + 6H{1 + r, 4, r, 0, x) + 3//(0, 1 + r, r, 0, x) 
-//(O, 0, 0, 0, x) - H{0, 0, 1, 0, x) + H{0, 1, 0, 0, x) - H{0, 1, 1, 0, x) 
+3H{1, 1+r, r, 0, x)-H{l, 0, 0, 0, x) -//(I, 0, 1, 0, x)+H{l, 1, 0, 0, x) 

-H(l, 1, 1, 0, x) - a/3 [(i/(r, 0, 0; 1) + 2H{r, 0; 1) + H{4, r, 0; l))i/(l, 

+//(r, 0; l)i/(0, 1, x) + i/(r, 0; 1)H{1, 1, x) 
2(4 - x) 



(387) 



a; 



xy^x(4 — x) 
+2H{4:,r,0,x) 



({2)H{r, x) - H{r, 0, 0, x) - H{0, r, 0, x) 



(388) 




2(4-D) /■ 
^ ' i=— 1 



12^13 



(389) 

(390) 



(l-x 



-i7(0,x), 



(391) 



C(2)ii-(1, x) + ii"(0, 1, 0, x) - 2i/(l, 0, x) + i/(l, 1, 0, x) 

[C(2) + 2i7(0, x) - H{0, 0, x) + H{1, 0, a;)] , (392) 
(l-x)L J 

^ ^amih ^) + 4//(l, 0, x) - 4C(2)//(0, 1, x) - 4C(2)i/(l, 1, x) 



a; 
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where: 



2 p(39) 







-2H{1, 0, 0, x) - H{0, 0, 1, 0, x) + H{0, 1, 0, 0, x) - 4H{0, 1, 1, 0, x) 
-H{1, 0, 1, 0, 1, 0, 0, x) - 4//(l, 1, 1, 0, x)^ 



(l-x) 



2C(2) 



+2C(3) + (4 - C(2)) i/(0, x) +C(2)i/ (1, - 2i/ (0, 0, x)+2H{l, 0, x) 
+i/(0, 0, 0, x) + H{0, 1, 0, x) - H(l, 0, 0, x) + H(l, 1, 0, x)] . (393) 




^ ' i=— 1 



1 

X 



H{0,x), 



(394) 
(395) 

(396) 



2^3 



(l-x) 
1 



+ 



H{r, 0; 1)//(1, x) + 3H{1, 0, x) - 2//(l + r, r, 0, x) 

C(2) + VSH{r, 0; 1) + 2H{0, x) - H{0, 0, x) + H{1, 0, x) 
3 



a/2^(4 - x) 



1 + 



(l-x). 



^(r,0,a;), 



(397) 



a; 



3C(2)iJ(l, x) + 6H{1, 0, x) - 3H{r, r, 0, x) - 2C(2)iy(l + r, r, x) 



+3H{0, 1, 0, x) - 3H{1, 0, 0, x) + 3i/(l, 1, 0, x) - 9H{1 + r, r, 0, x) 
+2/^(1 + r, r, 0, 0, x) + 2H{1 + r, 0, r, 0, x) - 4if(l + r, 4, r, 0, x) 

"2 



-2//(0, 1 + r, r, 0, x) - 2H{1, 1 + r,r,0,x) + Vs 



-H{r,0;l)H{0,l,x) 



+ '^H{r, 0, 0; 1)H{1, x) + 3H{r, 0; 1)H{1, x) + '^H{4, r, 0; 1)H{1, x) 



+ -H{r,0;l)H{l,l,x 
3 



2C(2) + 2C(3) 



(1 - x) 

+AH{0, x) - H{0, x)C{2) + H{1, x)Ci2) - 2H{0, 0, x) 
+2H{1, 0, x) + H{0, 0, 0, x) + H{Q, 1, 0, x) - H{1, 0, 0, x) 

+H{1, 1, 0, x) - 3H{1 + r, r, 0, x) + ^3 [i/(r, 0, 0; 1) + 2H{r, 0; 1) 
+H{r,0;l)H{l,x) + H{4,r,0;l) 



^/x{4: - x) 



(l-x) 



[C(2)if(r, x) + 2H{r, 0, x) - H{r, 0, 0, x) 



-//(0,r,0,x) + 2i/(4,r, 0,x)] . 



(398) 
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The following two diagrams are reducible to a combination of two different 4- 
denominator diagrams, as explained in Eqs. (5,6) of [l,^, by simple partial fractioning. 




^ J ViV^V, 



12-'-^13 



(399) 
(400) 



where: 




aF[f = 


1 

X 


aF[f = 


1 - 

X . 


aF[f = 


1 - 

X . 


aFr = 


1 - 

X . 



aFl 



(40) 



where: 

aF^ 



(401) 

l-i7(0,x)j, (402) 

1 - C(2) - i/(0, x) + if(0, 0, x) + if(l, 0, x)l , (403) 

1 - C(2) - 2C(3) - (l - C(2)) //(O, x) + C(2)i/(1, x) + i/(0, 0, x) 
-i7(0, 0, 0, x) + iJ(0, 1, 0, x) + if(l, 0, x) - iJ(l, 0, 0, x) 
+i/(l,l,0,x)], (404) 

1 - C(2) - ^C'(2) - 2C(3) - if(0, x) + (C(2) + 2C(3))//(0, x) 

+ (C(2) + 2C(3))//(1, x) + (l - C(2)) (if(0, 0, x) + H{1, 0, x)) 

+H{<d, 1, x)C(2) + C(2)i/(1, 1, x) - /J(0, 0, 0, x) + i/(0, 1, 0, x) 
-//(I, 0, 0, x) + H{1, 1, 0, x) + iJ(0, 0, 0, 0, x) + if(0, 0, 1, 0, x) 
-i/(0, 1, 0, 0, x) + if(0, 1, 1, 0, x) + H{1, 0, 0, 0, x) 

+if(l,0,l,0,x)-i/(l,l,0,0,x) + /7(l,l,l,0,x) . (405) 

(406) 
(407) 




^ ' i=—l 



H{0,x), 



aFr 



(41) 



(1-x) 

-if(l,0,x) + — ^ 
X [1 — x) 



C(2) + i7(0,x) 
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(408) 

i/(0,0,x) + i7(l,0,x) ,(409) 



aFl 



(41) 



2C(2)i7(l, x) + 2H{1, 0, x) + 2H{0, 1, 0, x) - 2H{1, 0, 0, x) 
1 



-2H{l,l,0,x) 



x) 



C(2) + 2C(3) + if(0,x)-C(2)/f(0,x) 



aK 



(41) 



+C(2)i^(l, x) - H{0, 0, x) + H{1, 0, a;) + H{0, 0, 0, a;) + if (0, 1, 0, x) 
-if(l,0,0,x) + if(l,l,0,x)] , (410) 

- [C(2)//(1, x) + 2C(3)if(l, x) + 0, x) - C(2)i/(1, 0, x) 
X L 

+C(2)/i(l, 1, x) + Ci2)H{0, 1, x) + H{0, 1, 0, x) - H{1, 0, 0, x) 

1, 0, x) + H{0, 0, 1, 0, x) - H{0, 1, 0, 0, x) + H{0, 1, 1, 0, x) 

+H{1, 0, 0, 0, x) + H{1, 0, 1, 0, x) - H{1, 1, 0, 0, x) + H{1, 1, 1, 0, x) 



(1-x) 



C(2) + ^e{2) + 2C(3) + (l - C(2) - 2C(3))i/(0, x) 

+ (C(2) + 2C(3))i/(l, x) - (l - C(2)) H{0, 0, x) + C(2)i^(0, 1, x) 

+H{1, 0, x) - C(2)i/(1, 0, x) + C{2)H{1, 1, x) + H{0, 0, 0, x) 
+iJ(0, 1, 0, x) - H{1, 0, 0, x) + H{1, 1, 0, x) - if(0, 0, 0, 0, x) 
+H{0, 0, 1, 0, x) - H{0, 1, 0, 0, x) + H{0, 1, 1, 0, x) + H{1, 0, 0, 0, x) 



+H{1, 0, 1, 0, x) - H{1, 1, 0, 0, x) + H{1, 1, 1, 0, x) 



(411) 
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